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SECTION  I 


INTRODUCTION 

1.  1  RELEVANCE  TO  AIR  FORCE  NEEDS 

Low  total  impulse  requirements  for  Air  Force  satellites  are 
presently  provided  by  small  monopropellant  rocket  engines  and  resisto- 
jets.  The  growth  in  total  impulse  capability  needed  for  the  near  future 
is  being  met  by  improving  the  useful  lifetime  of  these  engines.  However, 
future  satellites  requiring  much  larger  total  impulse,  either  for  attitude 
control  and  station  keeping  for  long  periods  or  for  larger  impulse  maneu¬ 
vers,  such  as  orbit  repositioning,  will  require  the  high  specific  impulse 
available  with  electric  propulsion  to  avoid  the  unacceptable  increases  in 
propellant  mass  of  the  lower  1^^  chemical  thrusters.  Development  of 
the  pulsed  ablative,  solid  propellant  thruster  and  the  colloid  thruster  is 
proceeding  for  the  relatively  low  impulse  stationkeeping  missions.  To 
perform  orbit  repositioning  with  lighter  weight  than  the  chemical  systems 
now  in  use,  and  to  provide  the  capability  for  future  requirements  of 
multiple  orbit  shift  maneuvers,  a  compact,  lightweight,  reliable  electric 
thruster  is  needed.  A  pulsed  electromagnetic  thruster  employing  a 
strong  bias  magnetic  field  is  an  attractive  candidate  for  meeting  this 
requirement  with  high  efficiency. 

1.  2  COMPARISON  OF  ELECTRIC  THRUSTERS 

The  performance  of  several  nonthermal  types  of  electric  thrusters 

capable  of  greater  than  1000  seconds  I  is  summarized  in  Figure  1. 

sp  j  2 

The  electrostatic  types  include  the  colloid  and  ion  thrusters.  '  The 

3 

electromagnetic  types  include  the  steady  flow  MPD  arc  thruster  and 

three  pulsed  types.  These  are  the  coaxial  thruster  using  injected  gas 

4  5  6 

propellant,  the  ablative  Teflon  thruster  '  and  the  pulsed  inductive 

7 

thruster.  The  latter  device  is  the  only  plasma  thruster  that  does  not 
use  electrodes;  the  pulsed  plasma  current  is  driven  inductively  and 
flows  in  a  closed  circular  path. 

It  is  interesting  that  the  performance  data  for  all  of  the  thrusters, 
except  the  colloid  engine,  group  fairly  well  around  a  straight  line  pass¬ 
ing  through  the  origin.  Furthermore,  except  for  the  low  value  of  the 
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Figure  1.  VVersas  i  for  Several  High  I 
®  sp  ®  sp 

Electric  Thrusters 

Convair  coaxial  device  and  the  high  values  of  one  of  the  Fairchild  Teflon 

thrusters,  all  of  the  plasma  thrusters  have  an  efficiency  at  any  l^p  that 

is  about  l/lOO  of  the  I  .  The  slope  of  n  versus  I  ^  is  a  little  more  than 
,  sp  sp 

twice  this  value  for  the  ion  engine  while  the  colloid  thruster,  at  1100 

seconds  I  ,  has  five  times  this  slope, 
sp 

The  practical  importance  of  the  ratio  H/l  is  that  it  is  proper- 

9  J-/ 

tional  to  the  thrust  per  unit  power.  If  m  is  the  rate  of  mass  flow  leaving 

2 

a  thruster,  then  l/2  m  (g  I  )  is  the  corresponding  rate  of  beam  kinetic 

sp 

energy,  where  the  effective  (thrust  producing)  velocity  is  g  1^^.  The 
total  power  into  the  thruster  is  therefore 


Power  = 


1/2  m  (g  I  )■ 


or,  since  the  thrust  is 


T  =  m  g  I 


Z 


T 

power 


and  the  thrust  available  from  a  fixed  power  source  increases  directly 

with  n/i  . 

sp 

For  stationkeeping  of  present  satellites,  the  thrust  available  from 
the  spacecraft  power  is  of  greater  importance  than  the  propellant  con¬ 
sumption  which  is  considerably  less  than  that  for  competitive  chemical 
thruster  stationkeeping  systems.  The  ion  engine  produces  about  half  the 
thrust  of  the  colloid  engine  from  the  same  amount  of  power  but  it  is 
capable  of  operating  at  both  higher  power  levels  and  higher  specific 
impulse  values.  It  is  therefore  a  possible  candidate  for  future  space¬ 
craft  missions  requiring  higher  power  than  is  now  available  and  larger 
total  impulse  than  the  present  stationkeeping  requirements.  Its  principal 
disadvantages  are;  (1)  possible  spacecraft  contamination  from  the  mer¬ 
cury,  or  cesium,  propellant,  (2)  large  size  relative  to  chemical  or 
plasma  thrusters,  and  (3)  complex  power  conditioning. 


The  motivation  for  the  present  plasma  acceleration  study  has  been 
to  demonstrate  the  feasibility  of  developing  large  thrust  levels  at  high  I 
with  a  ratio  of  v  Igp  comparable  to  the  ion  engine,  with  a  noncontaminat¬ 
ing  gaseous  propellant  and  with  a  relatively  simple  and  compact  power 
conditioner,  propellant  system  and  thruster.  An  attractive  system  con¬ 
cept  is  that  it  may  be  possible  to  operate  the  thruster  on  gaseous  com¬ 
bustion  products  obtained  from  an  on-board  chemical  rocket  propellant 
system.  Other  attractive  features  of  pulsed  plasma  propulsion  systems 
that  are  inherent  in  the  pulse  mode  are  the  thrust  variability,  accurate 
control  of  small  impulse  bits  and  accelerated  testing  of  components  at 
high  pulse  rates. 


The  practical  constraint  of  small  physical  size  and  weight  and 
low  average  power  level  limit  the  scope  of  accelerator  possibilities, 
especially  the  scale  of  the  devices.  The  TRW  thruster  indicated  in 
Figure  1  was  12  inches  in  diameter  and  this  was  thought  to  be  an  upper 
limit.  The  present  investigation  was  undertaken  in  an  effort  to  find  a 
more  efficient  acceleration  mechanism  than  the  pressure  of  a  magnetic 
field  confined  by  the  diamagnetic  effect  of  a  plasma  current  sheet. 


3 


It  was  hoped  that  acceleration  by  some  form  of  magnetic  field  annihila¬ 
tion  in  the  interaction  in  the  collision  of  two  current  sheets  might  serve 
this  purpose. 

This  report  describes  the  results  of  a  5-year  study  of  impulsive 
plasma  acceleration  with  several  devices  that  have  used  some  form  of 
magnetic  field  annihilation  to  produce  the  acceleration  force. 
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SECTION  II 
BACKGROUND 


The  accelerator  circuits  used  in  all  of  the  impulsive  plasma 
accelerators  investigated  in  this  program  have  been  inductively  coupled 
to  the  current  carried  by  the  plasma.  In  these  devices,  charge  stored 
at  high  voltage  in  an  energy  storage  capacitor  is  conducted  through  a 
coil  by  closing  a  spark  gap  switch.  During  the  discharge  period  of  a 
few  microseconds,  the  rapidly  changing  current  in  the  coil  induces  an 
oppositely  directed  current  in  the  plasma  which  is  physically  located  as 
close  as  possible  to  the  accelerator  coil.  The  repulsive  force  between 
the  two  currents  accelerates  the  plasma  away  from  the  accelerator  coil. 
When  two  coils  of  this  type  are  used,  they  can  be  arranged  so  that  the 
magnetic  fields  in  the  individually  accelerated  current  sheets  are  anti¬ 
parallel.  Consequently  the  magnetic  fields  diffusing  into  the  composite 
current  sheet  after  collision  are  annihilated,  giving  up  their  energy  to 
plasma  heating  and  acceleration. 

An  earlier  study  of  impulsive  acceleration  in  an  inductive  device 
had  been  done  using  a  single  flat  spiral  coil  that  made  one  or  two  turns 
around  the  coil  axis.  By  using  many  closely  spaced  wires  the  accelera¬ 
tor  current  approximated  a  uniform  disc  of  current  so  that  a  sheet  of 
nearly  constant  current  density  was  induced  in  the  plasma.  The  rising 
magnetic  field  surrounding  the  accelerator  coil  was  prevented  from 
expanding  away  from  that  coil  by  the  highly  diamagnetic  plasma  current 
sheet  so  that  the  resulting  magnetic  pressure  was  effective  in  accelera¬ 
ting  the  plasma.  The  acceleration  mechanism  in  this  type  of  accelera¬ 
tor  was  investigated  in  static  fill  operation  by  means  of  optical  diagnostics 
and  miniature  probes  which  measured  the  plasma  current  density  and 
electric  and  magnetic  fields.  It  was  shown  that  the  plasma  current  was 
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Turely  azimuthal  and  was  carried  predominantly  by  electrons.  The 
axial  electromagnetic  body  force  associated  with  this  current,  which 
provided  the  acceleration  impulse,  was  balanced  by  a  combination  of 
electron  pressure  gradient  and  an  axial  charge  separation  electric  field. 
Toward  the  front  of  the  current  sheet  the  negative  electron  pressure 
gradient  represented  a  positive  acceleration  force  acting  to  accelerate 


5 


the  electrons  ahead  of  the  current  sheet.  To  preserve  neutrality  of  the 
plasma,  a  charge  separation  electric  field  developed  which  transferred 
the  pressure  gradient  force  to  the  ions.  At  the  rear  of  the  current  sheet 
this  effect  was  reversed  but  deceleration  of  the  ions  was  prevented  by 
the  electromagnetic  body  force  on  the  current  carrying  electrons. 

Toward  the  center  of  the  current  sheet  the  pressure  gradient  force  was 
small  and  the  more  familiar  situation  was  encountered  where  the  body 
force  on  the  electrons  was  transferred  to  the  ions  by  the  charge  separa¬ 
tion  electric  field.  The  overall  axial  impulse  imparted  to  the  plasma 
was  that  produced  by  the  electromagnetic  force  alone  while  the  electron 
pressure  gradient  forces  served  only  to  distribute  the  impulse  through 
the  current  sheet  differently  from  the  body  force  distribution.  This 
detailed  knowledge  of  the  current  sheet  structure  influenced  the  selection 
of  the  two  flat  spiral  coils  for  the  first  field  annihilation  accelerator. 

Also  this  configuration  provided  good  accessibility  of  the  plasma  volume 
for  optical  and  probe  diagnostic  measurements. 
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SECTION  III 


SELF- FIELD  ACCELERATOR  EXPERIMENTS 

In  a  self-field  accelerator  the  magnetic  field  that  couples  with  the 
plasma  current  to  produce  acceleration  arises  solely  from  the  circuit 
and  plasm.a  currents.  These  fields  develop  rapidly  in  a  pulsed  accelera¬ 
tor,  providing  the  local  electric  fields  that  drive  the  plasma  current. 

Thus  this  current  is  related  to  the  rate  of  change  of  magnetic  field. 

Since  the  local  acceleration  force  is  the  product  of  this  rate  of  change 
with  the  total  magnetic  field  at  any  point,  the  local  body  force  is  related 
to  the  local  rate  of  change  of  magnetic  pressure.  This  had  not  been  found 
to  be  large  enough  for  efficient  acceleration  in  the  single  flat  spiral  coil 
accelerator  and  it  was  hoped  that  a  significant  improvement  would  result 
from  the  two  coil  experiments. 

A  wave  mechanism  for  conversion  of  magnetic  field  energy  to 
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plasma  kinetic  energy  was  described  theoretically  by  Petschek  and 
subsequent  experimental  work  by  Bratenahl  and  Yeates,  particularly 
schlieren  photographs  of  the  accelerated  plasma,  Figure  2,  appeared  to 
confirm  the  theoretical  model.  Bratenahl' s  accelerator  used  two  side- 
by-side  inverse  Z-pinch  discharges  to  produce  two  cylindrical  sheets  of 
outward  moving  plasma.  After  collision  of  these  sheets,  the  plasma  and 
field  geometry  resembled  that  depicted  in  Figure  3,  which  is  reproduced 
from  Reference  9.  The  neutral  line  at  y  =  0  in  Figure  3,  is  midway 
between  the  axes  of  the  inverse  pinches  For  the  two  flat  spiral  coil 
accelerator,  the  neutral  line  is  a  circle  at  about  80%  of  the  maximum 
coil  radius  and  the  region  outside  of  that  radius  also  approximates 
Petscheck's  field  configuration.  The  inner  region  is  a  shock  preheated 
S-pinch,  This  is  illustrated  in  Figures  4  and  5  which  show  the  accelera¬ 
ted  current  sheets  before  collision  and  the  field  annihilation  configuration 
that  developes  after  collision. 

The  second  magnetic  field  annihilation  self-field  accelerator  used 
two  flat  sheets  of  parallel  coil  wires  in  place  of  the  flat  spiral  coils. 

Both  of  these  self-field  accelerators  are  described  below. 


Figure  2 


Schlieren  Photographs  of  Argon  Shock  Collision  in 
Bratenahl's  Double  Inverse  Pinch  (Dark  circles  ar 
current  rods  spaced  10  cm  apart) 


Figure  3.  Petschek's  Magnetic  Field 
Annihilation  Configuration 


igure  4.  Magnetic  Fields  and  Current  Figure  5.  Magnetic  Fields  and  Current  Sheets 

Sheets  Before  Collision  During  Field  Annihilation 


.• 
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3.  1  FLAT  SPIRAL  COIL  FIELD  ANNIITILATION 
ACCELERATOR 

The  accelerator  eniployetl  two  40  cm  diameter  flat  spiral  coils  facing 
each  other  and  spaceti  closely  enough  that  each  coil  was  inductively  coupled 
to  the  plasma  current  at  any  location  in  the  space  between  the  coils.  Each 
coil  comprised  24,  one-turn  sprirals  of  No.  14  copper  wire  extending  from 
an  OD  of  40  cm  to  an  ID  of  10  cm.  The  wires  were  returned  radially  from 
the  inner  edge  of  the  spiral  to  the  outer  edge  where  connections  were  made 
to  2  meter  long  RG8-U  coaxial  cables  that  separated  the  coils  from  the 
capacitor  bank.  The  bank  consisted  of  six,  4.  5  pf  capacitors  with  a  sep¬ 
arate  spark  gap  connected  to  each  capacitor. 

Figure  6  is  a  front  view  of  one  of  the  coils  showing  the  spiral  wires 
on  the  front  surface  and  radial  wires  on  the  back  side  which  are  located 
2  inches  behind  the  spiral  wires.  The  coil  wires  were  covered  with  two 
10  mil  sheets  of  Mylar  and  a  1/16  inch  thick  annular  glass  disc  with  a 
7.  5  cm  diameter  hole  at  the  center.  The  axial  distance  between  the  two 
glass  covers  was  3.62  cm,  or  about  9%  of  the  coil  diameter. 

The  coils  were  mounted  in  a  small  vacuum  chamber.  They  were 
inserted  through  circular  holes  in  square  Lucite  side  plates  as  shown  in 
Figure  6.  These  plates  rested  against  an  open  steel  framework  \yith  solid 
upper  and  lower  steel  plates  and  glass  windows  supported  on  O-ring  seals 
on  the  sides.  A  side  view  of  the  assembly  is  shown  in  Figure  7.  The 
glass  covers  are  shown  extending  beyond  the  coil  structure  in  this  figure. 
They  were  later  replaced  by  covers  having  the  same  OD  as  the  coil. 

As  shown  in  Figure  7,  the  capacitors  and  spark  gaps  are  located 
below  the  coils.  A  closer  view.  Figure  8,  shows  the  cable  connections. 

The  white  tubes  shown  in  this  figure  circulated  dry  air  through  the  gaps. 

The  accelerator  was  operated  without  preionization  or  crowbarring 
A  static  fill  of  argon  gas  was  used  at  500  millitorr  pressure.  The  bank 
voltage  was  20  kv  (5400  joules).  At  this  voltage  and  pressure,  a  dia¬ 
magnetic  plasma  current  sheet  formed  with  no  observed  delay  as  can  be 
seen  from  the  circuit  current  (lower  trace)  in  Figure  9.  The  initially 
large  rate  of  current  rise  followed  by  a  reduced  rate  of  fall  at  the  end 
of  the  first  half-cycle  is  characteristic  of  close  coupled  inductive  accel¬ 
erators.  Another  characteristic  feature  of  the  circuit  current  in  pinch 
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Figure  8.  View  Showing  Cable  Connections  at 
Coil  and  Spark  Gap  Housings 

r’0vices  can  also  be  seen  in  this  figure.  At  about  2  psec  after  the  current 
starts  rising  in  each  half  cycle  there  is  a  kink  in  the  curve  where  dl/  dt 
suddenly  increases.  At  this  time  the  individual  planar  plasma  current 
sheets  have  collided  at  the  midplane  between  the  coils  and  the  impedance 
due  to  the  rate  of  inductance  increase,  which  was  positive,  suddenly  goes 
to  zero.  The  associated  voltage  drop,  I(dL/  dt)  must  then  be  supplied  by 
an  increase  in  the  other  inductive  voltage  drop  term,  L(dl/  dt).  The  upper 
trace  in  Figure  9  is  argon  II  light  viewed  along  the  axis  of  the  accelerator 
(three  shot  overlay).  A  photograph  of  the  accelerator  firing,  illuminated 
by  its  own  light,  is  shown  in  Figure  10. 

The  region  of  maximum  brightness  is  a  thin  annular  disc  where  the 
sheets  of  accelerated  plasma  collide  at  the  midplane.  This  is  seen  in 
the  time  exposure  photographs  shown  in  Figures  11  and  12.  This  zone 
appears  to  be  1  to  2  mm  thick  and  has  a  radial  extent  nearly  the  same  as 
that  covered  by  the  coil  wires. 
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Oscillograph  Traces  Showing  Circuit  Current 
(Tower)  and  Argon  II  Light  (Upper).  Three 
traci  s  overlay,  1  psec/  cm  sweep  speed. 

Figure  10.  Time  Exposure  of  Accelerator  Firing 


Figure  11.  Time  Exposure  of  Dis-  Figure  12.  Time  Exposure  of 

charge-Sidc  View  ,  Discharge-Inclined 

View  Showing  Dark 
Center  Disc 
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3.1.1  Magnetic  FiclH  Mapping 

Magnetic  field  menieasurements  were  made  with  miniature  probes 
located  at  the  tip  of  a  glass  tube  that  was  introducted  into  the  plasma  from 
the  top  of  the  chamber  and  extended  parallel  to  the  coil  surfaces.  The 
probe  shaft  therefore  ran  parallel  to  the  plane  of  the  plasma  current  sheet. 
The  probe  carriage,  which  permitted  continuous  positioning  both  axially 
and  vertically  (radial)  can  be  seen  at  the  top  of  the  vacuum  chamber  in 
Figure  6.  The  probe  itself  is  visible  in  Figures  1  1  and  1  2  where  it  is 
located  at  the  midplane  in  the  thin  region  of  maximum  brightness. 

A  check  on  possible  interference  from  the  probe  shaft  was  made  by 
two  sets  of  data,  one  with  a  6  mm  diameter  glass  tube  and  the  other  with 
a  3  mm  diameter  glass  tube.  The  probes  were  individually  calibrated 
and  the  measurements  were  made  two  months  apart.  Spot  checks  of  the 
results  showed  agreement  to  within  ±5%. 

The  probes  used  in  the  6  mm  diameter  glass  tubing  consisted  of 
20  turns  with  a  2  mm  diameter  window.  The  second  pair  of  probe  coils 
were  about  half  this  dimension  in  the  axial  and  circumferential  directions. 

In  order  to  retain  sensitivity  they  were  extended  about  8  mm  in  the  rf.dial 
direction.  This  was  felt  to  be  acceptable  since  radial  variations  in  this 
accelerator  geometry  are  small  compared  to  axial  ones  and  there  are  no 
gradients  in  the  circumferential  direction. 

The  coil  geometry  of  the  axial  magnetic  field  probe  was  a  rectangle 
1  mm  wide  in  the  circumferential  direction  and  9  mm  long  in  the  radial 
direction.  The  coil  had  10  turns  of  #36  wire.  The  radial  magnetic  field 
probe  was  a  one  layer  solenoid,  I  mm  in  diameter  by  8  mm  long  in  the 
radial  direction,  with  40  turns  of  #36  wire.  The  probes  were  calibrated 
with  Helmholtz  coils  driven  with  a  sinusoidal  generator,  the  current  being 
measured  by  a  Hewlett-Packard  Model  456  A  current  probe. 

Magnetic  probe  data  were  recorded  at  a  sweep  speed  of  0.5  psec/cm. 
The  probes  were  terminated  in  52  ohms  at  the  oscilloscope  and  their  sig¬ 
nals  integrated  passively  with  a  time  constant  of  100  psec.  Data  scans  were 
made  at  2  cm  radial  intervals  and  2  mm  axial  intervals.  The  resulting 
oscilloscope  traces  were  read  and  plotted  directly  as  functions  of  time 
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with  axial  position  as  parameter  for  each  radial  station.  Radial  magnetic 
field  plots  are  shown  in  Figure  13  through  16  for  radial  stations  of  8,  12, 

16  and  20  cm,  respectively. 

Axial  distributions  can  be  visualized  in  these  figures  as  the  inter¬ 
section  with  the  data  plots  of  vertical  planes  passing  through  the  constant 
time  lines.  For  times  up  to  about  1.  5  psec,  the  propagating  magnetic 
field  typical  of  the  single  coil  accelerator  can  be  seen.  About  1  psec  is 
required  for  the  toe  of  the  field  distribution  to  reach  the  midplane  at 
18.  1  mm.  At  about  1. 6  psec  the  field  lines  begin  to  dip  and  reach  a  min¬ 
imum  at  about  2  psec.  This  time  corresponds  to  the  sudden  rise  in  dl/dt, 
shown  as  the  kink  in  circuit  current  in  Figure  6,  that  occurs  when  plasma 
layers  collide  at  the  midplane. 

Another  interesting  feature  of  the  radial  magnetic  field  plots  is  the 
linear  decrease  in  with  z  after  the  separate  current  sheets  have  co¬ 
alesced  into  a  single  one.  Petschek's  illustration  in  Figure  3  shows  a 
finite  region  at  the  midplane  where  =  o  and  the  field  energy  stored  in 
has  already  been  converted  to  kinetic  energy.  The  present  data  do  not 
show  such  a  region.  decreases  linearly  with  z  all  the  way  to  the  mid¬ 
plane. 

The  general  magnitude  of  the  radial  magnetic  field  is  comparable  for 
radial  stations  of  8,  12,  and  16  cm,  but  has  decreased  noticeably  at  the 
20  cm  station.  The  16  cm  radial  station  is  the  neutral  point  (or  circular 
ring)  where  Bz  =  o.  It  is  not  only  the  neutral  point  of  the  fully  developed 
field  configuration  (usually  a  consequence  of  being  the  locus  of  initial 
collision  of  the  two  sheets)  but  is  also  the  location  of  Bz  =  o  for  the  separate 
axially  accelerating  current  sheets  as  a  consequence  of  the  coil  geometry 
which  produces  a  maximum  radial  magnetic  field  at  this  radius. 

The  magnetic  field  maps  in  Figures  13  through  16  provide  a  visualiza¬ 
tion  of  the  current  sheet  development  as  time  progresses.  Spatial  field 
distributions  can  be  made  for  any  time  by  erecting  lines  from  the  r-t  plane 
to  the  corresponding  contour  plots  and  drawing  a  continuous  curve  through 
the  intersection  points  .  Flux  plots,  or  maps  of  the  magnetic  field  lines, 
are  also  helpful  in  conveying  an  intuitive  feeling  of  the  accelerator  opera¬ 
tion.  For  example,  the  vacuum  field  in  the  space  between  the  coils  is 


15 


16 


essentially  axial,  so  that  deviations  from  an  axial  map  of  magnetic  flux 
contour  lines  (i.e,,  lines  of  constant  magnetic  flux)  are  seen  as  the  con¬ 
sequence  of  plasma  currents.  At  early  times,  the  plasma  current  sheets 
near  the  coil  confine  the  magnetic  field  to  the  layer  between  each  coil  and 
its  current  sheet.  Within  this  layer  the  resultant  field  is  nearly  radial 
and  the  thickness  of  the  layer  grows  with  time  as  the  current  sheet  moves 
toward  the  midplane.  Another  interesting  feature  is  the  x-type  neutral 
point  which  must  appear  at  a  radius  of  16  cm  and  must  have  the  same 
general  appearance  as  Figure  3,  The  development  of  this  x-type  field 
configuration  as  time  progresses,  as  well  as  the  progress  of  the  initial 
current  sheets  across  the  accelerator  and  the  diffusion  of  field  lines  ahead 
of  the  current  sheet  can  be  seen  in  field  line  maps. 

Such  maps  have  been  prepared  from  the  experimental  data  using  a 
time  sharing  computer  with  a  Cal  Comp  plotter  and  existing  subroutines 
that  accept  a  matrix  of  experimental  data  points,  read  from  hand  faired 
curves,  and  fits  curves  to  these  points  that  have  continuous  first  and  second 
derivatives.  This  made  it  feasible  to  compute  and  plot  current  density 
contour  maps  from  the  magnetic  field  data  in  addition  to  the  field  contour 
maps . 

The  spatial  intervals  of  the  original  probe  measurements  were  2  cm 
in  the  radial  direction  and  2  mm  in  the  axial  direction.  This  is  a  matrix 
of  1 1  0  test  shots  each  for  the  radial  and  axial  magnetic  fields.  These  220 
oscilloscope  traces  were  plotted  and  faired  to  make  consistent  families  at 
each  radial  station  like  those  shown  in  Figures  13  through  16.  These  data 
were  then  read  at  0.  1  psec  intervals  and  input  to  the  computer.  This  was 
a  total  of  40  time  steps  for  220  plots  or  8800  input  data  values.  The  final 
key  punched  values  were  then  printed  out,  examined  for  inconsistencies  and 
edited. 

The  three  characteristic  phases  of  this  accelerator  are:  first,  the 
period  of  axial  motion  when  each  current  sheet  is  essentially  unimpeded 
by  the  other  and  the  rising  inductance  dL/dt  describes  the  real  load  im¬ 
pedance  in  ohms,  second,  the  collision  phase  beginning  with  the  forward 
edge  of  the  current  sheets  which  carry  little  moving  plasma  and  ending 
with  the  collision  of  the  axially  moving  mass  layers,  and  third,  the  period 
of  radial  mass  motion  of  shock  heated,  highly  conductive,  plasma  which 
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carrit-’s  (lux  inward  inside  the  neutral  radius  and  expels  it  outward  at 
larger  radii.  These  jjliases  are  shown  in  the  three  groups  of  flux  contour 
maps  of  Figures  17a,  17b  and  17c,  respectively.  Figure  18  shows  a  series 
of  10  plots  of  both  current  density  and  magnetic  flux  contours  for  a  time 
span  of  4.0  psec  at  0.4  psec  intervals.  The  flux  contour  intervals  are 
0.  0005  Webers  for  all  frames  and  the  current  density  intervals  are  5x10 
amps/M^,  with  zero  cor r(*sponding  to  the  heavy  lines.  The  upper  plot  in 
each  frame  is  current  density  and  the  lower  plot  is  magnetic  flux. 


Figure  17a.  Magnetic  Flux  Lines 
from  0.  Z  to  1.2  psec, 
2  X  lO*!  W/Li  ne 


Figure  17b,  Magnetic  Flux  Lines 
from  1 , 4  to  2.  4  psec, 
3  X  10^  W/Line 


Figure  17c 

Magnetic  Flux  Lines  from 
2.  5  to  3.  0  psec,  3  x  10"^ 
W/  Line 
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18.  Plasma  Current  Density  and  Magnetic  Flux 
Contour  Maps  for  the  Time  Intervals 
Shown  (Continued) 


Figure  18.  Plasma  Current  Density  and  Magnetic  Flux 
Contour  Maps  for  the  Time  Intervals 
Shown  (Continued) 


Plots  of  current  density  versus  axial  position  for  16  cm  radius,  at 

several  times,  are  shown  in  Figure  19.  Except  for  the  4.0  psec  case, 

these  plots  correspond  to  the  contour  maps  shown  in  Figure  18.  At  1  .  2 

psec  the  sheets  have  begun  to  merge  at  the  midplane  and  at  2.  0  psec  the 
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current  density  has  risen  to  a  maximum  of  10,  000  amps/ cm  ,  or  about 
50%  more  than  the  maximum  of  the  individual  current  sheets.  At  this  time 
the  total  current  sheet  thickness  is  a  little  more  than  one  cm.  The  thick¬ 
ness  increases  slightly,  but  the  current  sheet  remains  concentrated  near 
the  midplane  as  it  decays.  It  is  interesting  that  the  current  density  max¬ 
imum  has  receded  from  the  midplane  at  2.  4  psec.  This  is  discussed  in 
more  detail  later. 


Figure  19.  Plasma  Current  Density  Distribution 
at  1 6  cm  Radius 


3.1.2  Plasma  Conductivity 

The  neutral  line  at  16  cm  radius  is  a  convenient  location  for  an 
experimental  determination  of  conductivity  because  the  mass  velocity  and 
both  magnetic  field  components  are  zero  at  all  times  so  that  the  conductivity 
is  simply 
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The  probe  measurements  already  described  determine  both  and 
;  the  latter  is  the  time  variation  of  flux  (B^  probe)  within  the  neutral 
radius.  However  rapid  changes  in  both  and  occur  during  collision 
so  that  better  time  resolution  is  required  for  a  good  ratio  determination 
than  was  possible  from  those  data.  Accordingly,  the  electric  field  was 
measured  with  a  probe  having  a  pair  of  exposed  electrodes  separated 

1  .  2  cm  in  the  0  direction.  Each  probe  tip  had  a  small  wire  (antenna) 
attached  to  it  to  increase  the  current  collection  area.  This  wire  extended 

2  mm  from  the  probe  tip  in  the  positive  and  negative  z  directions  (4  mm 
total  length).  This  introduced  some  uncertainty  in  the  E^  structure 
because  of  possible  averaging  over  the  center  4  mm  of  the  current  sheet. 

As  a  check  on  the  electric  field  probe  data,  Eq  was  also  determined  from 
magnetic  flux  variation  measured  with  the  probe  and  the  agreement 
was  reasonably  good. 

The  curl  of  the  magnetic  field  is  determined  almost  entirely  by  the 
radial  component  at  the  midplane.  The  spatial  variation  of  the  current 
density  can  therefore  be  accurately  determined  by  B^  probe  measurements. 
Figure  20  shows  the  current  density  distribution  over  the  central  4  mm  of 
the  current  sheet  for  times  ranging  from  1  .  2  to  3.  4  ps  in  0.  2  psec  steps. 

As  a  check  on  consistency,  the  average  current  density  over  the  Rogowsky 
window  area  was  determined  from  these  distributions  and  found  to  be 
satisfactory,  as  shown  in  Figure  21.  Also  shown  in  Figure  21  is  the  high 
resolution  current  density  at  the  midplane  (zso  in  Figure  19).  The  sudden 
drop  at  2.  3  psec  is  related  to  the  collision  of  the  plasma  masses.  This 
effect  also  appears  in  the  electric  field  which  is  compared  with  in 
Figure  22.  The  maximum  current  density  at  2,  2  psec  is  related  to  the 
maximum  electric  field  at  that  time. 

The  conductivity,  determined  from  Equation  (1)  is  also  shown  in 
Figure  22.  It  increases  about  50%  when  the  plasma  masses  collide  and 
remains  at  the  higher  level  as  radial  motion  develops. 
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Figure  20.  Axial  Distribution  of  Current  Density 
Near  Midplane 


t  (i-SfC) 


Figure  21.  Comparison  of  Current  Density  at 
Midplane  for  Rogowsky  and  Curl  H 
Measurements 


Figure  22.  Time  Variation  of  Conductivity  Deter¬ 
mined  by  Eg  and  Jg  at  Coil  Midplane 

3.1,3  Spectral  Measurements 

3 .  1 .  3 .  1  Electron  Temperature  at  Discharge  Center 

The  intensities  of  several  argon  11  {once  ionized  argon)  spectral 
lines  were  used  to  measure  the  electron  temperature  at  the  axis  of  the 
coils.  For  the  relatively  high  densities  typical  of  this  experiment,  the 
populations  of  the  higher  lying  atomic  levels  are  to  a  good  approximation 
controlled  by  electron  collisions  with  radiative  processes  being  relatively 
unimportant.^*  Since  electron  collisional  excitation  is  balanced  by  its 
inverse  process,  electron  collisional  deexcitation,  and  since  the  collision 
times  are  much  shorter  than  the  time  scale  of  the  experiment,  these 
atomic  levels  are  populated  according  to  a  Boltzmann  distribution.  Hence, 
the  intensity  of  line  radiation  from  these  levels  is  given  by 
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I  =  Chv'  g  A  e 


E 

KT 

e 


or 


In 


In  (Ch) 


(2) 


where 

I  =  line  intensity 
C  =  constant 
h  =  Planck's  constant 
V  ~  frequency  of  the  transition 
g  =  multiplicity  of  the  upper  state 
A  =  atomic  transition  probability 
E  =  upper  state  energy 
K  =  Boltzmann's  constant 


The  slope  of  the  line  determined  by  plotting  In 
experimental  values  for  I,  is  -  . 


— ; —  versus  E,  using 
g  A  V  ^ 


e 

The  intensities  of  seven  argon  II  lines  from  the  4S-4P,  3D-4P,  and 
4P-5S  arrays  were  measured.  These  lines  originate  from  levels  having 
energies  varying  from  19.2  (above  the  ion  ground  state)  to  24.6  ev  and 

O 

range  in  wavelength  from  about  4000  to  5000  A.  Light  from  the  plasma 
was  imaged  onto  the  input  slit  of  a  half  meter  Jarrel-Ash  monochromator 
by  a  200  mm  telephoto  lens.  The  line  of  sight  was  axial  and  passed 
through  the  center  of  the  discharge.  A  photomultiplier,  whose  output  was 
recorded  on  an  oscilloscope,  detected  light  at  the  exit  slit. 


The  result  of  the  measurement  for  t  =  5.  5  gsec,  the  time  of  peak 
continuum  emission  at  the  discharge  axis,  is  displayed  in  Figure  23. 
Transition  probabilities  listed  by  Wiese  et  al.  ^  ^  were  used  to  produce 
this  graph.  The  experimental  points  are  a  reasonable  fit  to  a  straight 
line  with  T^  equal  to  3.  8  ev. 
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UPPER  STATE  ENERGY  LEVEL  pv) 

Figure  23.  Spectral  Line  Intensity  Measurement 

3  1.3.2  Spectral  Line  Pro  file  Measurements  at  r  =  o 

The  measurement  of  an  argon  II  line  profile  at  the  discharge  center 
was  undertaken  to  determine  either  the  ion  temperature  or  electron  den¬ 
sity.  In  the  limit  of  low  electron  density  and  high  ion  temperature,  the 
principal  line  broadening  mechanism  is  Doppler  broadening  caused  by  the 
thermal  motion  of  the  emitting  ions;  if  the  ions  are  in  thermal  equilibrium, 
the  profile  has  a  Maxwellian  shape  whose  width  is  given  by  the  ion  tem¬ 
perature.  In  the  other  limit  where  the  electron  density  is  high  and  the 
ions  are  cold,  the  dominant  broadening  mechanism  is  Stark  broadening 
caused  by  the  perturbation  of  the  emitting  ions  by  the  fields  of  nearby 
particles.  The  Stark  broadened  line  shapes  of  heavy  elements  can  be 
closely  approximated  by  the  Lorentzian  function  ^  ^ 


The  half-width,  A\  ^  ,  depends  strongly  on  electron  density  and  very  little 
on  temperature.  ^ 


o 

Profile  measurements  of  the  argon  11  line  at  4806  A  were  made. 

The  experimental  setup  was  the  same  as  that  for  the  line  ratio  measure- 

O 

ment  at  the  axis  except  that  the  entrance  slit  was  adjusted  for  O.Z  A 


5.  4  jjLsec,  which  is  coincident  with  the  maximum  continuum  emission. 
Figure  25  shows  the  profile  at  5.4  psec  compared  with  a  Lorentzian  hav- 

O 

ing  a  full  width  of  1.  12  A  at  half  of  maximum  intensity.  The  excellent 
agreement  is  strong  evidence  for  Stark  broadening.  The  possibility  that 
the  line  profile  is  Doppler  broadened  but  that  the  ions  are  not  in  thermal 
equilibrium,  is  unlikely  but  cannot  be  entirely  ruled  out.  Assuming  an 
average  ion  energy  given  by  the  half  width  of  the  profile  and  the  electron 
density  equal  to  the  atomic  fill  density,  the  energy  equipartition  time  for 
ions  is  about  2  psec.  The  observed  magnetic  field  compression  along 
with  the  continuum  peak  implies  that  the  plasma  is  compressed  so  that  the 
energy  equipartition  time  is  probably  less  than  2  psec. 

Assuming  Stark  broadening  as  the  most  likely  mechanism,  it  is  pos¬ 
sible  to  infer  an  electron  density  from  the  theoretical  line  broadening 

parameters,  given  by  Griem,  and  the  measured  half  width  and  electron 

o  23  3 

temperature.  Taking  AA*  1  .  1  2  A  yields  n^*  4.  0x10  /M  atr»0  and 
t  a  5.4  psec.  An  electron  density  this  high  would  require  a  strong  com¬ 
pression  at  the  center  of  the  discharge  or  a  more  modest  compression  and 


Figure  25.  Argon  II  Spectral  Line  Profile  at  5.4  psec 
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multiple  ionization.  Since  multiple  ionization  would  probably  produce  A 

O 

111  line  emission,  a  search  for  the  A  III  line  at  3858  A  was  undertaken; 

however,  line  emission  at  3858  A  could  not  be  detected.  For  single  ioniza- 
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tion  and  n  *4x10  /M  ,  it  would  be  necessary  to  compress  the  plasma 
e 

initially  in  the  region  inside  of  the  neutral  radius  into  a  disk  about  2  cm 
wide  by  3  cm  in  radius.  Framing  pictures  showed  a  luminous  region  with 
similar  dimensions  but  with  a  1  cm  radius  hole  or  dark  region  at  the 
center. 

To  reiterate,  the  observed  argon  II  line  broadening  is  most  likely 

caused  by  Stark  broadening.  The  density  implied  by  the  half  width  is 
23  3 

4  X  10  /M  .  Doppler  broadening  from  non- Maxwellian  ions  seems 
unlikely  because  of  the  short  energy  equipartition  times. 

3.  1 . 3.  3  Spectral  Line  Profiles  at  the  Midplane 

In  an  attempt  to  observe  directed  radial  motion  of  the  plasma,  by 
means  of  Doppler  shifts  of  spectral  lines,  line  profiles  were  measured 
along  a  radial  line  of  sight  in  the  midplane  of  the  discharge.  The  line  of 
sight  through  the  center  was  blocked  at  r  =  0  so  that  light  from  only  the 
half  diameter  towards  the  observer  was  accepted.  If  this  were  not  done, 
outward  radial  motion  for  example,  would  be  observed  as  both  a  blue 
shift  of  light  emitted  by  ions  closest  the  observer  and  a  red  shift  of  light 
emitted  by  ions  on  the  far  side  of  the  plasma.  The  unavoidable  perturba¬ 
tion  of  the  obstructing  object  made  the  data  for  the  time  of  the  compres¬ 
sion  at  the  center  not  valid. 

The  experimental  procedure  was  similar  to  that  described  for  the 
line  profile  measurements  at  r  =  0.  In  addition,  the  center  wavelengths 
of  the  lines  measured  at  the  midplane  were  compared  with  the  wavelengths 
of  the  same  lines  originating  near  one  coil  face  at  the  beginning  of  the 
discharge.  Such  line  radiation  emanating  from  near  the  coils  during  the 
initial  gas  breakdown  should  not  be  Doppler  shifted  because  plasma  motion 
does  not  have  time  to  occur.  The  positions  of  the  line  centers  of  the  two 
sets  of  lines  agreed  within  instrumental  precision  or  about  0.  1  A. 

O 

Figure  26  presents  profiles  for  the  argon  II  line  at  4848  A  while 

o 

Figure  27  shows  the  profiles  for  the  argon  III  line  at  3858  A.  A  salient 
feature  of  these  spectra  is  the  drop  in  line  emission  at  the  time  the  cur¬ 
rent  sheets  collide.  The  continuum  radiation  shows  a  burst  at  this  time. 
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Figure  26.  Argon  II  Spectral  Line  Profiles 

One  possible  explanation  is  that  when  the  current  sheets  collide  the  ions 
are  excited  to  the  higher  levels  of  ionization  (i.  e.  ,  triple  ionization  or 
higher)  causing  the  argon  II  and  III  line  emission  to  be  quenched.  It  would 
be  simple  to  test  this  hypothesis  by  searching  for  argon  IV  line  emission; 
however,  no  argon  IV  lines  are  available  in  the  visible  spectrum.  The 
argon  in  line  peaks  at  2.  5  psec  while  the  argon  II  maximum  emission 
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Figure  27.  Argon  III  Spectral  Line  Profiles 
occurs  at  3  to  4  jisec.  The  plasma  is  probably  twice  ionized  at  2.  5  psec 

e 

since  the  3858  A  argon  III  line  emanates  from  a  level  29.  7  ev  above  the 
argon  III  ground  state  or  56.  1  ev  above  the  argon  II  ground  state.  For 
an  electron  temperature  of  about  4  ev  it  is  unlikely  that  this  level  would 
be  populated  enough  for  observable  emission  if  the  ions  were  singly  ionized 


.  0 

Figure  28  shows  the  profile  of  the  4848  A  argon  II  line  at  3  psec 

0 

along  with  a  Lorentzian  with  a  half  width  of  0,  34  A  for  comparison.  The 

close  agreement  between  the  experimental  points  and  the  comparison 

curve  is  strong  evidence  for  Stark  broadening.  If  the  broadening  para- 

meters  given  by  Griem  are  used  the  inferred  electron  density  is  1. 2  x 
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10  /M^.  This  is  consistent  with  a  compression  of  the  initial  neutral 

fill  density  into  an  8  mm  wide  band  at  the  midplane  with  twice  ionized  ions. 
The  framing  pictures  show  that  the  luminous  region  of  the  discharge  (the 
luminosity  of  heavy  ion  plasmas  is  usaully  associated  with  electron  den¬ 
sity),  is  confined  to  such  a  band. 


Figure  28.  Argon  II  Spectral  Line  Profile  at  3  psec 

3.1.  3.4  Electron  Temperature  Determination  at  the  Neutral  Radius 

Measurement  of  the  electron  temperature  at  the  neutral  radius 

allows  a  comparison  of  the  measured  conductivity  with  that  inferred  from 

1  5 

the  classical  (Spitzer)  relation.  Anomalously  low  conductivity  may 
result  from  micro-instabilities  having  electric  fields  which  scatter 
electrons  and  produce  a  large  effective  electron  collision  rate.  Since 
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nonclassical  conductivity  has  been  observed  in  the  earlier  field  annihila¬ 
tion  experiment,  it  is  important  to  compare  the  conductivity  in  this 
experiment  with  that  inferred  from  classical  theory. 

Again,  the  electron  temperature  was  determined  from  measured 
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argon  II  line  intensities.  The  high  electron  diensities  (about  10  '  /M^) 
typical  of  this  plasma  at  the  midplane  after  the  sheets  collide  insure  that 
partial  thermodynamic  equilibrium  holds  for  the  atomic  transitions  used. 
That  is,  collisional  processes  as  opposed  to  radiative  processes  deter¬ 
mine  the  populations  of  the  upper  states  of  the  optical  transitions  mea¬ 
sured.  Five  argon  II  lines  were  measured  along  a  line  of  sight  in  the 
midplane  of  the  discharge  passing  through  a  minimum  radius  equal  to 
9  cm.  This  line  of  sight  sampled  light  emitted  from  plasma  with  an 
average  radial  position  equal  to  that  of  the  neutral  radius.  Poor  optical 
access  through  the  coils  precluded  a  measurement  in  the  axial  direction 
which  would  have  given  better  radial  resolution.  Time  did  not  permit 
measurement  at  several  radii  which  would  be  needed  for  an  Abel  inversion 
to  obtain  radial  resolution. 

The  monochromator  entrance  slit  width  was  set  for  an  instrumental 
resolution  of  1.  3  A.  This  was  adequate  since  the  most  broadened  lines 

O 

had  widths  which  were  typically  0.  4  A.  If  the  slit  width  had  been  made 
wider  to  include  the  far  wings  of  each  line,  the  increased  continuum  radi¬ 
ation,  falling  within  the  instrument  band  width  would  have  made  the  line 
intensity  measurement  more  difficult.  As  it  was  the  strong  continuum 
level  made  it  impossible  to  measure  the  line  intensities  accurately  enough 
to  infer  electron  temperature  for  the  tjmes  of  peak  continuum  emission, 
i.  e.,  from  2  to  2.  7  psec. 

Figure  29  shows  the  results  of  the  temperature  measurement  as  a 
function  of  time.  The  temperature  at  1, 5  psec  is  about  3.  7  ev,  slowly 
increasing  to  approximately  5  ev  at  4  psec. 

Figure  30  shows  the  calculated  value  of  or,  using  T^  and  the  Spitzer 
relation,  compared  with  the  value  measured  from  E_  and  j  .  The  ion 

23  3  D  b 

density  has  been  taken  to  be  .10  /M  .  From  1  to  2  psec  the  agreement 
is  excellent,  while  at  3  psec  the  calculated  value  is  about  25%  less  than 
the  measured  value.  No  temperature  measurement  was  possible  for  the 
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T  e  vs  T  IME 
AT  NEUTRAL  RADIUS 


Figure  29.  Electron  Temperature  at  the  Neutral 
Radius 
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TIME  (MSEC) 

Figure  30.  Conductivity  at  the  Neutral  Radius 
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interval  from  Z  to  2.  7  ^sec.  It  is  important  to  note  that  the  measured 
conductivity  after  collision  is  larger  than  the  Spitzer  value.  If  there  were 
an  additional  resistive  mechanism,  the  measured  value  would  be  smaller, 
not  larger,  than  the  classical  value.  Thus  the  comparison  demonstrates 
that,  in  this  experiment,  the  conductivity  is  determined  solely  by  clas¬ 
sical  Coulomb  collisions. 

3.1.4  Time  Resolved  Photography 

3 .  1.4.  1  Center  Region 

Figure  31  shows  image  converter  framing  camera  pictures  at  the 
discharge  center.  The  line  of  sight  is  axial  and  the  region  seen  has  its 
maximum  radius  equal  to  4.  1  cm. 

For  heavy  ion  plasmas,  such  as  this  one,  the  plasma  luminosity  is 
closely  identified  with  plasma  density.  At  3.  5  psec  the  plasma  has  moved 
into  the  central  region  and  has  an  inner  radius  of  2.  5  cm.  At  5.  5  psec  the 
plasma  is  fully  compressed  into  an  annular  region  bounded  by  0.6  and  3  cm 
radii.  Magnetic  probe  measurements  show  that  flux  which  leaks  through 
the  poorly  conducting  gas  at  the  center  at  early  times  is  strongly  compres¬ 
sed  by  the  inward  moving  plasma,  seen  in  these  photographs.  The  exclu¬ 
sion  of  the  moving  plasma  by  this  compressed  field  produces  the  hole  seen 
at  the  center.  Subsequent  to  this  time  the  plasma  decays  away  as  can  be 
seen  in  the  pictures  for  6.  5  and  7.  5  psec. 

Figure  32  shows  a  streak  photograph  of  the  center  of  the  discharge 
using  a  mask  aligned  horizontally  and  the  line  of  sight  parallel  to  the  axis 
of  the  device.  The  time  axis  is  up,  and  the  streak  duration  is  from  0  to 
10  psec.  The  implosion  can  clearly  be  seen.  The  inward  velocity  of  the 
plasma  front  is  inferred  to  be  about  1  cm/ psec. 
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f  -  3.5  seconds 
t  =  4.5  seconds 
t  =  5.5  /i seconds 
t  =  6. 5  p seconds 
^  -  7.5  n  seconds 

Figure  31  .  Framing  Pictures  of 

the  Discharge  Center  Figure  32.  Streak  Picture  of  the 
Viewed  Axially  Discharge  Center 

from  t  =  0  to  10  p-sec 

Figure  33  shows  a  streak  for  the  interval  1.  5  to  3,  5  psec  after  the 
discharge  is  initiated.  At  about  2  psec  a  burst  of  light  (seen  on  the  origi¬ 
nal  Polaroid  photograph)  is  emitted  by  the  plasma  which  persists  for  about 
0.2  psec.  This  is  at  the  same  time  and  duration  as  the  intense  continuum 
emitted  as  a  result  of  the  current  sheet  collision.  There  are  two  possible 
explanations.  First,  a  small  fraction  of  the  continuum  is  scattered  by 


Figure  33.  Streak  Picture  of  the 
Discharge  Center 
t  »  1 .  5  to  3.  5  psec 


8.5  m 
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solid  surfaces  along  the  line  of  sight  of  the  camera.  This  seems  unlikely 
because  both  surfaces  arc  far  removed  from  the  plasma  volume.  Second, 
the  plasma  at  the  center  is  excited  or  photoionized  by  the  UV  component  of 
the  continuum  and  subsequently  emits  line  radiation  as  it  is  deexcited. 
Strongly  accelerated  argon  plasmas  have  been  observed  to  be  excellent  UV 
sources  for  photoionizing  argon.  For  example,  30%  photoionization  of 
neutral  argon  by  an  inverse  pinch  discharge  was  observed  by  one  of  the 
authors  in  an  carlii  r  study.  It  should  be  noted  that  Rayleigh  or  Thompson 
scattering  of  the  continuum  by  the  plasma  at  the  center  (probably  not  fully 
ionized)  is  highly  unlikely  because  of  the  small  Thompson  and  Rayleigh 
scattering  cross  sections  (o-y  ~  10  while  ~  10  cm  ). 

3 .  1 . 4.  Z  Edge  View 

Figure  34  shows  the  center  region  of  the  discharge  viewed  along  a 
radial  line  of  sight.  The  center  is  necessarily  viewed  through  the  outer 
region  of  plasma.  At  1.25  psec  the  two  current  sheets  have  lifted  off  the 
coils  and  are  approaching  the  midplane  while  at  2.  25  psec  they  have  fully 
coalesced  and  the  luminosity  is  confined  to  a  region  about  4.  5  mm  wide. 

At  3.  25  psec  the  current  layer  has  broadened,  and  its  visible  output 
has  diminished  considerably.  The  radial  compression  has  become  apparent 
by  4.  25  psec  while  5.  25  psec  is  close  to  the  time  of  peak  compression. 

The  luminous  area  near  the  coils  is  the  second  current  sheet  which  forms 
as  the  voltage  at  the  coils  changes  sign. 


Figure  34 

Framing  Pictures  of  the 
Discharge  Center 
Viewed  Radially 
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An  estimate  of  radial  impulse  was  made  by  using  the  midplane 
measurements  of  and  B  and  assuming  that  the  corresponding  force 
density  acted  over  a  1  cm  wide  region  around  the  midplane. 


Figure  35  shows  radial  distributions  of  at  different  times.  The 
corresponding  current  density  map  is  shown  in  Figure  36.  The  resulting 
force  density  distribution  is  plotted  in  Figure  37.  The  total  radial  force 
for  the  inner  and  outer  regions  is  plotted  in  Figure  38  as  a  function  of 
time  and  the  corresponding  impulse  values  are  shown. 


R  Icm  I 


Figure  35.  Axial  Magnetic  Field  at  Coil  Midplane 
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Figure  37.  Radial  Force  Density  Distribution 
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Figure  38.  Time  Variation  of  Radial  Force 


The  efficiencies  and  mean  plasma  velocities  were  low.  The  highest 

4 

velocity  was  0.7  x  10  M/sec  for  the  inner  region.  This  was  about  half 
of  the  axial  velocity  before  collision.  For  the  outer  region  the  velocity 

4 

was  0.55  X  10  M/sec.  The  corresponding  efficiencies  were  2.1  and  1.4% 
respectively. 

3.1.6  Plasma  Collision  at  Neutral  Radius 


The  neutral  radius,  where  the  axial  magnetic  field  is  zero,  is  about 
16  cm.  The  radial  variation  of  this  magnetic  field  component  is  plotted 
for  different  times  in  Figure  39.  For  very  early  times,  the  neutral  radius 
is  near  the  outer  edge  of  the  coil  (20  cm).  It  moves  rapidly  inward  to 
16  cm  at  2.0  psec  and  then  moves  out  to  nearly  17  cm  after  2,6  psec. 


Figure  40  shows  the  variation  of  the  radial  magnetic  field  with 
axial  position  at  different  times  for  the  same  radius.  At  1 . 8  and  2.  0 
psec  the  curves  show  a  structure  which  seems  to  correspond  to  the  trailing 
edge  of  the  mass  layer.  The  slope  at  the  midplane  determines  the  current, 
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Figure  39.  Radial  Variation  of  Axial  Magnetic  Field  at  16  cm 
Radius  for  Different  Times 
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Figure  40.  Axial  Variation  of  Radial  Magnetic  Field  at  16  cm 


Radius  for  Different  Times 


which  is  plotted  as  a  function  of  time  in  Figure  41.  The  contribution  from 
the  other  component  of  the  curl  has  been  omitted  since  it  is  less  than  3% 
at  all  times. 

Figure  41  also  shows  the  time  variation  of  the  azimuthal  electric 
field  and  the  conductivity  determined  from  the  ratio  of  the  current  density 
to  the  electric  field. 


Oscilloscope  traces  of  the  radial  magnetic  field  are  shown  in 
Figures  42  and  43  for  axial  positions  at  0.  5  mm  intervals  near  the  midplane. 

All  of  the  traces  in  Figures  42  and  43  have  the  same  characteristic 
feature.  The  rising  magnetic  field  seen  at  any  fixed  axial  position  reaches 
a  maximum  after  which  it  rapidly  decreases.  A  few  tenths  of  a  micro¬ 
second  later  the  field  abruptly  stops  falling  and  remains  nearly  constant. 
Careful  examination  of  the  traces  shows  the  first  abrupt  change  to  be 
propagating  inward  toward  the  midplane  while  the  second  one  propagates 
outward  away  from  the  midplane. 

A  cross  plot  of  the  inward  moving  feature  shows  its  velocity  to  be 

4 

1,3  X  10  M/sec,  which  is  a  reasonable  value  for  the  axially  moving 
plasma  before  collision.  Thus  the  rising  field  before  the  maximum  is 
reached,  is  due  to  motion  of  the  radial  magnetic  field  which  accumulates 
within  the  mass  layer  because  of  the  high  conductivity  (the  diamagnetic 
effect)  and  the  subsequent  rapid  decay  corresponds  to  arrival  of  the 
region  behind  the  mass  layer  where  the  curl  of  the  magnetic  field  is  very 
low.  The  fact  that  the  peak  is  sharply  defined  implies  a  clean  separation 
between  the  trailing  edge  of  the  mass  layer  and  the  low  current  region 
behind  it. 


The  outward  moving  feature  propagates  more  rapidly  than  the 
earlier  inward  motion  of  the  mass  layer.  For  example,  the  end  of  the 
period  of  rapidly  decaying  field  for  the  lower  trace  in  Figure  43  occurs 
only  0.  1  fisec  before  the  corresponding  event  for  the  upper  trace  in 
Figure  42.  This  feature  is  believed  to  be  a  magnetoacoustic  wave  whose 
speed  is  given  by  the  relation 


a2 


+  yk  T  (n.  +  n  ) 
e  1  e 


n.M. 

X  1 


(5) 


45 


Variation  of  Current  Density,  Azimuthal  r-ie 
and  Conductivity  with  Time  at  16  cm  Radius 


Radial  Magnetic  Field  vs  Time  at  16  cm  Radius 
{z  =  2,  2,5,  3,  3.5,  4  and  4.  5  mm  from  midplane, 
0.5  dsec/cm  sweep,  lower  trace  is  circuit  current) 


Figure  43.  Radial  Magnetic  Fielfl  vs  Time  at  16  cm  Radius 
( z  =  0.  5,  1  ,  1.5  mm,  0.  5  [iscc/  cm  sweep, 
lower  trace  is  circuit  current) 

where  Y  is  the  specific  heat  ratio  5/3  and  k  is  Boltzmann's  constant. 

A  cross  plot  of  the  outward  moving  feature  shown  in  Figure  44  gives 

a  velocity  variation  with  axial  position  which,  with  the  corresponding 

magnetic  field  values  from  Figure  40  giv'es  a  constant  ion  density  of 
20  3 

about  0.5  X  10  per  M  for  all  positions  more  than  0.5  mm  away  from 
the  midplane . 
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Figure  44.  Arrival  Time  for  Outgoing  Waves  (R  =  16  cm) 
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Thus  nearly  all  of  the  mass  corresponding  to  the  initial  fill  density 
22 

of  1.6  X  10  is  compressed  within  a  1  mm  thick  disc  at  the  midplane  at 
an  average  density  of 


n. 

1 


3 . 6  cm 
0.  1  cm 


X  1. 6  X  10^^ 


5.8  X  10^^  per 


A  rough  check  on  electron  temperature  and  wave  speed  can  be 
made  from  the  curve  drawn  in  Figure  44.  The  minimum  slope  of  the 
arrival  curve  is  the  thermal  acoustic  speed  given  by  Equation  (5)  with 

4 

=  0.  The  value  shown  is  0.5  x  10  M/sec,  or  0.5  mm  per  0.  1  psec. 
This  is  consistent  with  an  ionization  level  of  2.  8  and  T^  =2.  8  ev. 

The  detailed  features  shown  in  Figure  41  for  both  Eg  and  ,  and 
consequently  a,  were  accurately  determined  by  the  oscilloscope  traces. 
The  current  density  was  determined  by  cross-plotting  the  traces  shown 
in  Figures  42  and  43  and  Eq  was  recorded  directly  as  illustrated 
in  Figures  45  and  46  for  two  radii. 

As  an  aid  in  interpreting  the  events  at  the  midplane  which  are 
shown  in  Figure  41,  the  conductivity  and  the  plasma  state  relation  (ion 
density,  plasma  temperature  and  degree  of  ionization  Z)  were  calculated 
for  local  thermodynamic  equilibrium  and  using  the  following  recent  data 
for  incremental  ionization  energies: 


Vj  =  15.75  ev 
-  27. 62  e V 
=  40.90 ev 
=  59. 79  e V 

The  result  is  plotted  in  Figure  47.  Also  shown  in  Figure  47  is  a  grid  of 
recombination  rates  calculated  according  to  Spitzer's  relation  for  radi¬ 
ative  recombination. 

The  most  striking  feature  of  Figure  41  is  the  abrupt  drop  in  conduc¬ 
tivity  associated  with  the  collision  of  the  incoming  mass  layers  at  2.  1  psec. 
The  0.  1  psec  delay  between  maximum  and  Eq  maximum  is  real  and 
must  be  accounted  for  in  an  explanation  of  the  event.  It  should  also  be 
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Figure  45.  Time  Variation  of  Eq  at  Midplane  (Upper  Trace) 
and  Circuit  Current  for  16  cm  Radius  -  0.  5  psec 
Sweep 


Figure  46.  Time  Variation  of  Eq  at  Midplane  (Upper  Trace) 
and  Circuit  Current  for  17  cm  Radius  -  0.  5  psec 
Sweep 
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Figure  47.  Equilibrium  Properties,  Radiative  Recombination 
Rate  and  Conductivity  for  Argon 


noted  that  the  outward  traveling  magnetoacoustic  wave  originates  at 
2,  2  psec  when  jg  reaches  a  minimum  and  the  break  in  the  downward 
slope  of  Eq  occurs.  It  is  also  noted  that  there  are  two  events  (1.4  and 
2.  1  psec)  producing  abrupt  loss  in  conductivity  and  each  is  accompanied 
by  an  abrupt  increase  in  the  slope  of  Eq. 

The  highest  conductivity  occurs  at  1.4  psec,  just  before  the  first 

4 

sudden  drop.  The  value  at  this  time  is  a  =  1.95  x  10  per  f2M.  According 
to  Figure  47  this  requires  both  high  temperature  and  high  density. 

The  spectroscopic  determination  of  temperature  indicated  values  as 
high  as  5  ev  at  later  times.  Values  at  least  this  high  must  also  exist  near 
1.4  psec  to  explain  the  measured  conductivity.  If  the  midplane  density  is 
assumed  to  increase  from  the  fill  density  (1.6  x  10^^)  to  10^^,  the  temper¬ 
ature  would  be  5.  5  ev  and  the  plasma  would  be  3.  5  times  ionized.  This 
density  increase  represents  a  six-fold  compression. 
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The  events  at  1 . 4  and  2.  1  [j-sec  are  regarded  as  the  arrival  of  two 
different  mass  layers  corresponding  to  doubly  and  singly  charged  ions. 
The  ratio  of  arrival  times  is  about  \2  as  would  be  expected  according  to 
this  argument. 

After  collision  the  conductivity  can  be  decreased  by  cooling  and  by 
rarefaction  (density  reduction  due  to  bounce-back  after  collision).  It  is 
interesting  to  examine  these  possibilities  for  the  two  collisions. 

It  is  assumed  on  the  basis  of  the  relative  losses  of  conductivity  in 
Figure  41  that  the  mass  involved  in  the  second  collision  is  twice  that 
involved  in  the  first  one.  A  six-fold  density  increase  at  1 , 4  usee 
corresponds  to  compression  of  the  fill  gas  within  6  mm  of  the  midplane 
to  a  final  thickness  of  1  mm.  Since  the  coil-to-coil  spacing  was  36  mm, 
this  means  that  1/6  of  the  initial  mass  has  arrived  to  within  1  mm  of  the 
midplane  at  1.4  psec.  At  this  time  1/3  of  the  remaining  mass,  car¬ 
ried  by  the  fast  mass  layers,  reaches  the  midplane  and  the  other  2/3 
arrives  at  2.1  psec.  In  terms  of  the  total  mass  initially  between  the 
coils  at  16  cm  radius,  17%  is  compressed  to  the  midplane  by  1.4  Usee 
onds,  28%  more  collides  with  this  17%  at  1.4  psec  and  the  remaining 
55%  arrives  at  2.  1  Usee . 

From  Spitzer's  equilibration  times,  ion-ion  equilibration  occurs  in 
less  than  1  nanosecond  and  ion-electron  equilibration  in  less  than  10  nano¬ 
seconds.  Consequently  equilibrium  can  be  assumed  for  all  of  the  data. 

The  equilibrium  properties  just  after  collision  at  1.4  psec  are 
calculated  from  an  energy  balance  that  equates  the  thermal  and  ionization 
energies  of  both  the  initial  and  incoming  masses,  plus  the  kinetic  energy 
of  the  incoming  mass,  to  the  thermal  and  ionization  energies  of  the  mass 
after  collision.  The  results  in  Table  1  were  obtained  from  this  calculation; 

Table  1.  Plasma  Properties  at  Midplane  for  16  cm  Radius  and  1.4  psec 


"i 

T 

Z 

a 

(per  M^) 

lev) 

(per  DM) 

Incoming  mass 

10^^ 

3.0 

2.0 

1.3  X  lo'^ 

Before  collision 

10^3 

5.5 

3.5 

1.95  X  10** 

After  collision 

1.5  X  10^^ 

5.2 

3.5 

1.95  X  lO"^ 

After  rarefaction 

4  x  10^^ 

4.6 

3.0 

1.5  X  10^ 
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The  ion  density  is  found  as  a  dependent  variable  in  this  calculation. 

The  energy  conservation  approach  defines  a  relation  between  plasma 

temperature  and  ionization  level  after  collision.  This  relation,  with  the 

one  specified  by  Figure  47,  determines  unique  values  of  Z  and  T  for  the 

4 

observed  conductivity  of  1.5  x  10  .  The  ion  density  is  then  read  from 
the  location  on  the  map  determined  by  Z  and  T. 

The  ion  density  after  collision  is  25%  of  that  before  collision  in 
spite  of  the  additional  mass  carried  by  the  colliding  layers.  This  is  con¬ 
sistent  with  the  picture  of  a  strong  rarefaction,  due  to  bounce  back  after 
collision,  which  launches  the  magnetoacoustic  wave  detected  by  the  mag¬ 
netic  field  probe. 

The  reduction  in  midplane  temperature  shown  in  Table  1  occurs 
without  energy  loss  (i.e.,  without  radiative  cooling).  It  results  princi¬ 
pally  from  the  increased  ionization  level  of  the  incoming  mass. 

The  same  analysis  was  made  at  2. 1  psec  with  the  results  shown  in 
Table  2. 

For  this  case  the  ion  density  is  very  low  after  collision  indicating 
a  much  stronger  rarefaction,  which  is  consistent  with  the  appearance  of 
a  stronger  outgoing  Alfven  wave  at  this  time  than  at  1.4  psec.  It  is  this 
rapid  decrease  in  ion  density  caused  by  the  bounce-back  after  collision 
that  accounts  for  the  large  and  sudden  loss  in  conductivity. 

The  current  density  maximum  in  Figure  41  is  reached  just  before 
the  rarefaction  begins.  The  sharpness  of  the  peak  indicates  a  sharp 


Table  2.  Plasma  Properties  at  Midplane  for  16  cm  Radius 
and  2.  1  psec 


n. 

1 

T 

Z 

0 

(per  M^) 

(ev) 

(per  QM) 

Incoming  mass 

10^^ 

3.0 

2.0 

1.3  X  lo'^ 

Before  collision 

2.6  X  10^^ 

3,5 

2.3 

1. 6  X  10^ 

After  collision 

1.5  X  10^^ 

3.6 

2.5 

1.5  X  10“* 

After  rarefaction 

4.5  X  10^^ 

3.0 

2.6 

0.8  X  10^ 
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density  gradient  at  the  trailing  edge  of  the  mass  layer.  The  reason  for 
this  is  that  the  mass  arriving  from  1. 95  to  2.  05  [isec  supplies  the  momen¬ 
tum  influx  needed  to  maintain  high  density  at  the  midplane  during  collision. 
When  this  momentum  influx  ceases  the  rarefaction  develops  radidly.  The 
maximum  in  Eq  occurs  later  than  the  Jq  minimum  because  of  increased 
magnetic  field  diffusion  speed  at  the  lower  conductivity  in  the  rarefied 
plasma. 

After  the  rarefaction  is  over,  conditions  at  the  midplane  from  2.  2 
to  2.  7  psec  resemble  those  during  the  initial  compression  from  0.  8  to 
1.4  psec.  The  conductivity  again  rises  resulting  in  an  increase  in  cur¬ 
rent  density  even  though  E0  continues  to  decrease.  The  plasma  is  ohm- 
ically  heated  to  about  3.  8  ev  according  to  the  spectroscopic  data  and  the 

rising  conductivity  indicates  one  and  one-half  decades  increase  in  ion 

23  3 

density  to  about  3x10  M  at  2.  7  psec.  As  the  ion  density  increases 
the  electron  density  e-fold  rate  rises  to  about 

1  ^”e  6 

-  —  — TT-  =  10  per  second 
n  dt  ^ 

e 


The  refore 


dn  / 

e  1  nO 

“JT*  =  -  n  10 
dt  e 


and  assuming  an  average  radiated  energy  of  30  ev  per  recombination, 
the  power  radiated  in  the  free-bound  continuum  is 


P  =  30  (1.6  X  10’^^)  10^  n  =  4.8  X  lO"^^  n 

e  e 
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or  for  n^  =  3  X  10  M  at  maximum  compression, 


P  =  1.44  X  10^^  watts/M^ 


This  is  four  times  the  ohmic  power  at  this  time. 
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The  radiating  volume  at  this  time  is  a  1  mm  thick  annulus  with  a 

-  5  3 

15  cm  ID  and  a  35  cm  OD,  or  8  x  10  M  .  The  radiant  power  is  there¬ 
fore  113  X  10  watts  at  2.7  psec  when  the  recompressed  plasma  reaches 
maximum  conductivity. 

According  to  Table  1  and  Figure  47,  the  e-fold  rate  at  the  first 
conductivity  peak  (1.4  psec)  is  also  about  10  indicating  a  comparable 
continuum  radiation  at  that  time.  However,  when  the  more  massive  cur¬ 
rent  sheet  arrives,  the  estimated  conditions  shown  in  Table  2  indicate 
a  recombination  rate  less  than  half  the  values  at  1 . 4  and  2.7  psec.  Yet 
it  is  at  this  time  (close  to  2  psec)  when  the  intense  continuum  radiation 
was  observed  in  the  spectroscopic  and  streak  data. 

It  follows  that  the  ion  density  immediately  after  collision  must  be 
about  a  decade  larger  than  the  value  shown  in  Table  2,  with  an  e-fold 
rate  of  about  3x10°  per  second.  This  suggests  that  the  bounce-back 
following  the  fast  ion  collision  at  1 . 4  psec  returns  in  phase  with  the 
slower  ion  collision  at  2.  1  psec  to  produce  a  higher  density  than  was 
estimated  by  the  steady  inflow  model  used  to  calculate  the  conditions 
listed  in  Table  2. 

3.1.7  Radial  Plasma  Motion 

Radial  plasma  motion,  inferred  from  arrival  times,  was  detected 
by  both  the  azimuthal  and  axial  electric  field  probes.  The  Eq  feature  con¬ 
sisted  of  a  propagating  step  and  the  latter  was  a  sharply  defined  spike  in 
E  as  illustrated  in  Figure  48.  Arrival  times  for  both  features,  plotted 
in  Figure  49,  indicate  that  they  belong  to  the  same  radially  propagating 
structure.  This  structure  is  defined  at  a  radius  as  small  as  17.5  cm 
(only  1 .  5  cm  outside  of  the  neutral  radius)  where  it  already  has  a  high 
initial  velocity.  This  velocity  is  constant  over  a  2  cm  wide  annular  re¬ 
gion  (from  16.5  to  18.5  cm  radius)  and  for  0.  3  psec  (from  1.6  to  1.9  psec). 

The  possibility  of  this  event  being  an  Alfven  wave  is  ruled  out  by  a 
consideration  of  the  axial  magnetic  fields  shown  in  Figure  39.  For  a  con¬ 
stant  wave  speed  from  17.9  to  18.6  cm  radius,  the  ion  density  would  be 
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required  to  increase  from  4x10  to  1  .  5  x  1  0  per  M  .  Both  the  magni¬ 
tude  and  the  radial  variation  are  far  from  plausible.  Evidently  this  radially 
propagating  feature  represents  a  mass  motion  of  substantial  velocity  that 
develops  before  collision  of  the  slow  more  massive  layers. 

An  estimate  of  the  radial  density  distribution  can  be  obtained  from 
the  relation 


je 

*^6  e  E 

r 

The  probe  data  required  for  this  calculation  are  plotted  in  Figure  50  at 
20  cm  radius  and  the  resulting  electron  density  distribution  is  shown  in 
Figure  51 . 

The  times  for  maximum  and  minimum  electron  density  are  plotted  in 
Figure  49.  The  latter  time  agrees  well  with  the  arrival  of  maximum 
as  should  be  expected  since  n^  is  a  minimum  when  E^  (in  the  denominator) 
is  a  maximum. 

The  possibility  that  the  maximum  and  minimum  events  in  Figure  51 
correspond  to  collision  and  bounce-back,  as  observed  at  16  cm  radius, 
is  ruled  out  by  the  fact  that  collision  of  the  massive  layers  (based  on 
extrapolation  of  collision  times  at  other  radii)  should  occur  at  2.6  llsec 

i.e,,  near  n^  minimum. 

It  seems  fairly  certain  that  a  thin  high  density  ring  of  plasma  is 

4 

mow.ng  radially  outward  at  3.2  x  10  M/sec  in  the  outer  few  centimeters 
of  the  accelerator.  It  is  about  3.3  mm  thick  (in  the  radij.1  direction).  If 
this  ring  were  3.3  mm  thick  in  the  axial  direction  also  it  would  contain 
about  1/4  of  all  the  initial  fill  mass  in  the  region  between  16  and  20  cm. 

The  kinetic  energy  of  this  mass  would  be  56  joules,  or  4%  of  the  capa¬ 
citor  bank  energy. 

There  are  two  acceleration  mechanisms  that  might  account  for  this 
effect.  One  is  the  radial  pressure  gradient  of  the  heated  ions  and  the 
other  is  the  Jq  electromagnetic  force.  An  order  of  magnitude 


56 


estimate  can  be  made  to  determine  whether  either  of  these  mechanisms 
is  large  enough  to  account  for  the  observed  radial  velocity. 

For  the  electromagnetic  force  to  produce  the  observed  velocity  at 

a  radius  as  small  as  17.5  cm,  the  acceleration  would  have  to  occur 

between  16  and  17.5  cm  where  B  has  the  required  direction.  However, 

z 

is  nearly  zero  in  that  region  for  the  first  lisecond  at  which  time  the 

plasma  has  advanced  about  1  cm  axially.  Therefore  only  0.4  psec  is 

available  for  outward  acceleration  while  the  plasma  moves  from  Z  =  1  cm 

to  the  midplane  at  Z  =  1.8  cm  at  a  speed  of  1.8  cm/psec.  An  upper 

limit  of  0.04  for  B  is  assumed  based  on  data  at  18  cm,  and  a  correspond- 

z  g 

ing  upper  limit  for  jo  of  10  is  used  based  on  data  at  16  cm.  Assuming 

°  22  3 

the  ion  density  to  be  equal  to  the  fill  density,  n.  =  1.6x10  /M  ,  the 

radial  acceleration  cannot  be  greater  than 


je  _  ' _ 0.04x10^ 

n.  M.  ~  7i  I  ,^22.  ,  _  , 

1  1  (1. 6  X  10  )  6,7  X  10 


3.7  X  10*^  M/sec^ 


4 

At  the  end  of  0.4  psec,  this  produces  a  velocity  of  only  0.  15  x  10  M/sec, 
or  about  20  times  less  than  the  observed  velocity.  Consequently  the  elec¬ 
tromagnetic  contribution  to  radial  acceleration  is  ignorable  for  a  time  as 
early  as  1.6  psec. 

An  upper  limit  can  be  placed  on  thermal  pressure  gradient  accel¬ 
eration  by  using  values  from  Table  1  just  before  collision  where  n.  = 

23  3  23  3  ^ 

1  0  /M  and  the  electron  density  is  3,  5  x  1  0  /M  ,  giving  a  total  parti- 

23  3 

density  of  4.  5  x  10  /M  .  This  electron  density  compares  well  with 
value  calculated  from  the  probe  data  at  20  cm  radius,  shown  in  Fig- 
51  . 

An  estimate  of  the  density  gradient  can  be  made  by  assuming  it  to 

i^^qual  to  that  corresponding  to  the  forward  slope  of  the  distrubution  in 

are  51.  The  leading  edge  passes  a  fixed  position  (i.e.,  the  density 

23 

creases  from  0.4  to  4,  0  x  10  )  in  0.  1  psec  at  a  velocity  (from  Fig- 

4 

re  49)  of  3.  2  X  10  M/sec.  The  corresponding  thickness  is  therefore 
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6  =  Vt  =  3.2  X  (10"^)  =  0.32  mm 
The  pressure  gradient  is  therefore 


dr 


kT  ^  =  (1.38  X  10"^^)(5.5  x  11605) 


23 

4.5  X  10 
0.32  X  10"^ 


1.24  X  10^ 


The  radial  acceleration  is 


dr 

^  ^  n.  M. 
1  1 


1.24  X  10^ 

- ZT - 

10  M. 


1.85  X  10^^ 


4  -7 

giving  a  velocity  , of  1.85  x  10  after  10  seconds. 

This  is  about  60%  of  the  observed  radial  velocity,  indicating  that 
the  radial  plasma  expulsion  is  due  to  thermal  pressure  stresses  and  not 
to  electromagnetic  forces. 

The  appearance  of  the  pressure  gradient  force  as  the  important 
accleration  mechanism  for  early  times  in  this  magnetic  field  annihilation 
experiment  was  unexpected.  Ordinarily  electromagnetic  forces  in 
plasmas  are  much  greater  than  pressure  forces.  In  this  device,  how¬ 
ever,  two  features  combine  to  produce  an  acceleration  similar  to  a 
shaped  charge  explosion.  One  is  the  high  density  produced  in  the  axial 
pinch  near  the  coil  midplane  and  the  other  is  the  progressively  slower 
axial  speed  of  the  plasma  outside  of  the  neutral  radius.  The  first 
feature  produces  high  density  and  the  second  produces  a  high  radial  den¬ 
sity  gradient. 

An  independent  measurement  of  the  plasma  radial  velocity  was 

O 

made  by  observing  the  Doppler  shifted  wavelength  of  the  4806  A  line. 

The  results  are  shown  in  Figure  52  and  the  arrival  time  and  velocity  are 
plotted  in  Figure  49.  The  agreement  with  the  other  data  is  satisfactory. 


59 


(V 


2 


V 

(i^) 


0  12  3 

t  (>iSEC) 


Figure  52,  Radial  Velocity  at  22  cm  Measured  by 
Argon  4806  Doppler  Shift 

An  autocollimator  was  used  for  this  measurement  with  the  line  of 

sight  tilted  20  degrees  from  the  radial  direction  and  passing  through  the 

/ 

midplane  between  20  and  22  cm  radii.  Initial  attempts  using  radial  line 
of  sight  gave  skewed  distributions  without  a  shift  in  the  wave  length  for 
maximum  intensity.  This  was  because  of  the  strong  background  of  un¬ 
shifted  light  within  the  accelerator.  It  was  eliminated  by  using  a  par¬ 
tially  masked  aperture  to  exclude  all  light  near  the  midplane  within  the 
accelerator  when  the  optical  axis  was  tilted  20  degrees. 

_*7 

The  initial  fill  mass  between  16  and  20  cm  radius  is  4.4  x  10"  kg. 

It  was  estimated  that  25%  of  this  mass  is  expelled  by  the  early  pressure 
gradient  force.  Assuming  all  of  the  mass  is  expelled,  the  remainder  is 

3 

ejected  by  the  subsequent  electromagnetic  force  at  9.  1  x  10  N-sec  im¬ 
pulse  (Figure  38).  The  results  listed  in  Table  3  are  then  determined. 

The  pressure  gradient  effect  contributes  almost  half  as  much  impulse  as 
the  electromagnetic  force  in  the  outer  region. 

The  detailed  study  of  the  current  sheet  collision  and  plasma  accelera- 
ation  in  the  two  flat  spiral  coil  accelerator  has  provided  a  good  understand¬ 
ing  of  its  operation.  From  an  initial  neutral  fill  density  of  1 . 6  x  10^^/M^, 

23  3  ’ 

ion  densities  of  the  order  of  3  x  1  0  /M  were  developed  first  in  a  very 

thin  annular  sheet  by  the  collision  and  later  near  the  axis  by  the  subsequent 

0-pinch.  Plasma  conductivities  ranged  from  to  2  x  per  with 
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Tabu-  3.  Comparison  of  Radial  Mass  Expulsion  by  Thermal  and 
Electromagnetic  Eorces 


Force 

Mass 

(kg) 

Mean 

Velocity 

(M/sec) 

Lnpulse 

(N-seconds) 

Kinetic 

Energy 

(Joules) 

%  Bank 
Energy 

The  rmal 

1.  1  x  10" 

7 

3.2  X  10^ 

3.5  X  10"^ 

56 

1.0 

Electro¬ 

magnetic 

3.3  X  10" 

7 

2.7  X  10^ 

9  X  10"^ 

1  20 

2.  1 

Both 

4.4  X  10" 

7 

2.  8  X  lo"^ 

12.  5  X  10"^ 

172 

3.  0 

corresponding  electron  temperatures  ranging  from  about  3  ev  before  col¬ 
lision  to  5  ev  after  collision.  The  collision  launched  a  pair  of  Alfven 
waves  that  traveled  from  the  midplane  back  toward  the  coils  and  in  the 
region  outside  of  the  neutral  radius  at  16  cm,  a  powerful  thermal  (pres¬ 
sure  gradient)  stress  produced  an  explosive  outward  acceleration.  These 
phenomena  can  be  understood  in  terms  of  the  current  sheets  commonly 
encountered  in  this  type  of  pulsed  inductive  accelerator.  At  least  in  the 
present  accelerator  where  the  acceleration  distances  were  of  the  order  of 
three  times  the  current  sheet  thicknesses,  the  wave-type  of  magnetic  field 
annihilation  acceleration  was  not  observed.  For  much  larger  ratios  of  these 
lengths  (i.e.,  magnetic  Reynold's  number),  as  occur  in  solar  flarec  the 
acceleration  processes  might  be  qualitatively  different.  For  the  present 
purpose  it  seemed  advisable  to  continue  the  investigation  with  an  accel¬ 
erator  designed  to  enhance  the  magnetic  acceleration  observed  in  the  outer 
region. 

3.  2  TWO-DIMENSIONAL  CURRENT  ACCELERATOR 

The  next  experiment  used  a  device  that  could  be  altered  in  various 
ways  to  enhance  the  magnetic  field  component  orthogonal  to  the  plasma 
current.  Three  coils  wound  on  lucite  blocks,  shown  in  Figure  53,  were 
used.  The  positions  of  all  of  the  coils  could  be  varied.  For  the  arrange¬ 
ment  shown  in  Figure  53,  the  large  coils  contribute  little  to  the  self-field 
required  for  acceleration.  This  field  component  could,  however,  be 
enhanced  by  tilting  the  large  coils  outward  with  the  small  coil  placed  at 
the  narrow  end,  as  shown  in  Figure  54.  The  reason  for  tilting  the  coils 
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Figure  53.  Coil  Configuration  for  Field  Annihilation 
Accelerator  (Current  density  probe  is 
shown  in  plasnia  region) 


Figure  54.  Schematic  of  Wedge  Shaped  Accelerator  Showing  Three 
Separate  Coils,  Plasma  and  Magnetic  Field 
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is  to  enhance  the  orthogonal  component  of  the  magnetic  field  so  that  the 
current  sustained  by  the  annihilation  of  the  antiparallel  axial  field  as  it 
diffuses  into  the  plasma  can  produce  a  greater  electromagnetic  force. 

The  coils  were  enclosed  by  a  sheet  of  copper  wrapped  closely 
around  them,  as  shown  in  Figure  55,  to  confine  the  magnetic  fields  with¬ 
in  the  plasma  volume.  This  is  illustrated  in  the  cross-sectional  view  of 
the  complete  accelerator  shown  in  Figure  56.  The  current  in  each  of  the 
large  coils  induces  a  countercurrent  that  flows  in  the  copper  sheet  and 
through  the  plasma.  Thus,  in  spite  of  the  fact  that  the  current  passing 
through  the  plasma  terminates  on  metal  surfaces,  this  device  is  an 
impulsive  inductive  accelerator.  Current  is  switched  into  the  coils  by 
a  spark  gap,  and  the  current  passing  through  the  plasma  runs  in  closed 
paths  around  the  coils.  These  paths  comprise  the  plasma  plus  the  metal 
sheets  wrapped  closely  around  the  outside  of  the  device. 

As  in  all  such  devices,  the  plasma  current  sheets  form  near  the 
coil  surfaces  and  move  away  under  the  action  of  the  electromagnetic 
forces.  The  time  rate  of  inductance  increase  related  to  this  motion  is 
the  useful  load  impedance  in  ohms  of  the  accelerator  circuit. 


Figure  55.  Photograph  of  Field  Annihilation  Accelerator  During 
Firing  (Three-Coil  Configuration) 


PLASMA 

CURRENT 


COIL 

CURRENT 


RETURN  CURRENT 
IN  METAL 

Figure  56.  View  Looking  into  Accelerator  Exit 
Showing  Currents  in  Coil,  Plasma 
and  Metal  Enclosure 

The  orthogonal  magnetic  field  component  has  been  calculated  and  is 

shown  in  Figures  57  through  59  for  a  range  of  coil  separations  and  tilt 

angles.  The  field  is  expressed  in  dimens’onless  form  by  using  the  axial 

length  of  the  accelerator,  X  ,  and  circuit  current,  I  .  This  is  the  total 

rn  C 

current  running  along  any  one  of  the  walls  and  it  is  the  same  for  each  of 
the  walls. 

The  upper  and  lower  coils  in  these  figures  have  six  16-gauge  wires 
connected  in  parallel  making  a  total  of  nine  turns  each  over  the  axial 
length  of  the  coil.  The  unshielded  inductance  of  each  coil  is  6  pHenrys 
When  the  outer  surface  of  the  coil  is  closely  covered  by  a  sheet  of  metal, 
the  inductance  is  a  function  of  x,  the  distance  between  the  coil  and  the 
plasma,  as  shown  in  Figure  60.  For  a  current  sheet  velocity  away  from 
the  surface  of  1.5  cm/psec,  the  motional  impedance  caused  by  the 
motion  for  a  single  coil  is 

^  =  =  (0.75  X  10-^)  J.5x  10^ 

=  1.1  ohms 

where  L  is  the  inductance  of  the  coil  when  the  plate  (which  simulates  the 
plasma  current  sheet)  is  at  the  distance  x. 
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I 


Figure  59.  Orthogonal  Magnetic  Field  on  Axis  of  Field 

Annihilation  Accelerator  (y  /X  =  0,20) 

'  o  m 


Figure  60.  Inductance  vs  Simulated  Plasma  Distance  for 
Large  Rectangular  Coil 
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Two  4.  5  fiF  capacitors  at  10  kv  were  connected  to  each  of  the  coils. 
For  a  mean  circuit  inductance  of  about  2  pHenrys,  this  gives  a  maximum 
circuit  current  of 


I 

m 


'9  X  10 


2  X  10 


A  A  v  A&x'.  v> 


Since  the  coil  has  nine  turns,  1  =  9  1  =  190,  000  amps.  For  a  unit 

cm  ^ 

value  of  the  dimensionless  field  shown  in  Figures  57  through  59,  i.  e.  , 


t:X  B 
m 


1. 


1.0 


and  for  a  15 -cm  long  coil  (x^  =  0.  15  M) 


B 


y  iTx 


(1.0) 


m 


.-7 


.  O.BWeber/M^ 


which  is  comparable  to  the  vacuum  ma  letic  fields  encountered  in  the 
flat  spiral  coil  accelerators. 

Quantitative  experiments  were  carried  out  with  the  two  large  coils 
tilted  as  shown  in  Figure  61,  The  accelerator  is  shown  during  operation 
in  Figure  62,  This  figure  also  shows  one  of  the  probes  that  was  mounted 
above  the  accelerator  on  a  traversing  mechanism  that  could  be  moved 
along  the  three  coordinate  axes.  Lateral  traverses  showed  uniform 
conditions  over  about  80%  of  the  length  of  the  coils.  Consequently,  all 
measurements  were  made  midway  between  the  side  plates,  with  Y  and  Z 
in  Figure  61  as  the  only  space  variables. 
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Conventional  miniature  probes  were  used  to  measure  the  plasma 
current  density  and  magnetic  fields.  The  Rogowsky  probe  (mentioned  in 
Figure  21  and  shown  in  Figure  52)  provided  a  direct  reading  of  current 
density  averaged  over  the  middiameter  of  the  toroidal  probe  coil.  An 
overall  check  on  probe  calibration  and  the  Rogowsky  technique  was  made 
for  this  accelerator  by  independently  determining  the  current  density 
from  the  curl-H  relationship  which,  because  of  the  relatively  weak  or¬ 
thogonal  fields,  depended  only  on  the  slope  of  the  antiparallel  Z  field  com¬ 
ponents  in  the  Y  direction,  i.e.. 


1 

■^x  ~  F  dY 

As  can  be  seen  in  Figure  63,  the  agreement  is  good  for  the  first  5  psec, 
after  which  the  probe  reading  is  too  low,  as  is  also  shown  in  Figure  21  . 
Thus,  impulse  in  the  Z  direction,  computed  from  measured  by  the 
Rogowsky  probe  and  B  from  a  magnetic  field  probe,  is  a  little  conserva- 
tive. 


Figure  63.  Comparison  of  from  Rogowsky  and  Magnetic 
Field  Probes  (z  =  6  cm,  y  =  0) 

Figure  64  shows  the  first  step  in  computing  efficiency  from  the 
probe  data.  The  current  and  field  probe  measurements  for  different  Z 
positions  are  multiplied  together  and  plotted  against  time.  Each  of  the 
resulting  curves  is  equivalent  to  an  oscilloscope  trace  that  would  be 
obtained  with  a  force  density  probe  if  such  a  thing  existed.  The  symbols 
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Figure  64.  Force  Density  vs  Time  at  Different  /ixial 
Stations  (Midplane) 


identify  the  Z  positions  as  well  as  the  times  for  points  read  from  the 
oscilloscope  traces  of  and  The  area  under  each  curve  is  the 

impulse  density  for  that  Z  position.  These  values  are  then  plotted 
against  Z  to  give  the  distribution  of  impulse  density  over  the  length  of 
the  accelerator  as  shown  in  Figure  65.  The  length  of  the  coil  is  also 
indicated  in  the  figure.  The  acceleration  force  is  negative  up  to  Z  =  5  cm 
with  a  negative  meiximum  at  2  cm.  The  positive  meudmum  occurs  at 
about  13  cm  and  the  useful  force  extends  6  cm  beyond  the  end  of  the 
accelerator. 


Figure  65  represents  an  upper  limit  on  impulse  in  that  the  probe 
measurements  along  the  Z  axis  (Figure  64)  were  assumed  to  act  over  the 
entire  cross-sectional  area  at  each  Z  location.  Even  with  this  assumption 
the  efficiency  is  low. 
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Figure  65.  Axial  Impulse  Density  Distribution 
The  efficiency  is  written  again  as 

1  2  , 

^  CV  2  ^ 

c  o 

where  m  is  propellant  mass,  v  is  plasma  velocity,  C  is  capacitance, 
is  initial  capacitor  voltage,  A1  is  the  impulse  bit  (mv)  and  is  the 
energy  initially  stored  in  the  capacitor.  For  this  experiment  the  capaci¬ 
tance  was  18  uF  and  V  was  10  kv.  The  total  mass  of  argon  in  the 
volume  between  the  coils  was  1.  i  x  lO"  kg  (corresponding  to  0.  5  torr  of 

_  3 

argon).  Thus,  the  efficiency  based  on  the  net  impulse  of  4.  92  x  10 
N-seconds  (Figure  65)  was 

n.  (4.92XI0-V 

2(900)(1.  1  X  10  ) 
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at  an  I 


sp 


of 


I  =  =  447  seconds 

sp  mg 

If  all  of  the  impulses  were  included  instead  of  the  net  value  only,  then 

n  =  ^  — W  =  5.0% 

2(900)(1.  1  X  10  ) 

at  an  I  of  900  seconds.  Even  at  this  ratio  of  n  to  1  ,  the  efficiency  pro- 

sp  sp 

jected  at  the  same  slope  to  1500  seconds  for  example,  would  be  only 

8.3%  as  compared  to  the  average  of  15%  shown  in  Figure  1. 

Probe  traverses  were  also  run  in  the  Y  direction  at  Z  stations  of 
2  and  12  cm  to  determine  the  deviation  between  the  Y- averaged  impulse 
density  and  the  midplane  value  at  y  =  0.  The  average  value  was  54% 
for  the  2  cm  station  and  87%  for  the  12  cm  position.  Thus,  the  above 
conclusions  based  on  the  assumption  of  constant  force  density  over  the 
cross  section  are  valid. 

The  experimental  studies  of  fairly  diamagnetic  plasmas  accelerated 

by  the  pressure  of  an  expanding  magnetic  field  have  shown  that  high  cur- 

2 

rent  densities  (of  the  order  of  10,000  amps /cm  )  and  good  conductivity 

4 

(better  than  10  per  ilM)  are  readily  achieved.  In  the  field  annihilation 
geometry  even  better  performance  is  seen  as  far  as  plasma  current  sheet 
properties  are  concerned.  But  the  conversion  of  the  field  energy  to 
kinetic  energy  is  consistently  poor  for  both  accelerators. 

The  fundamental  difficulty  with  the  field  annihilation  accelerator 
appears  to  arise  from  the  high  conductivity  of  the  current  sheet.  This 
can  be  seen  by  referring  to  Figure  66  which  shows  the  time  history  of 
the  plasma  current  density  and  the  magnetic  field  components  at  Z  =  1  2 
cm  on  the  midplane,  where  relatively  strong  acceleration  was  observed. 

At  the  midplane,  B  =  0  by  symmetry  and  the  rapid  buildup  of  B 

z  z 

in  the  Y  direction  seen  in  the  lov'er  trace  reflects  the  diaunagnetic  prop¬ 
erty  of  the  current  sheet.  It  is  the  upper  two  traces  in  Figure  66  (j^  and 
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cm 


Comparison  of  Orthogonal  Magnetic  Field  B, 
Plasma  Current  Density  at  x  =  y  =  0  and' 
=  12  cm  During  Current  Pulse.  is  also 
shown  for  0.  5  cm  increments  from  y  =  0  to 
y  =  3.5  cm  at  x  =  0  and  z  =  12  cm. 
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By)  that  show  why  high  conductivity  results  in  poor  efficiency  for  the  field 
annihilation  accelerator  (the  two  traces  are  an  overlay  of  two  shots  at 
different  oscilloscope  gains).  A  local  peak  in  both  the  orthogonal  field 
and  the  current  density  occurs  at  6  psec  and  the  corresponding  acceler¬ 
ation  is  a  maximum.  The  vacuum  field  corresponding  to  this  acceler¬ 
ator  geometry  only  arrives  after  7  psec  when  the  plasma  current  density 
is  rapidly  decaying.  After  4  to  5  psec  a  strong  plasma  current  has  devel¬ 
oped,  but  the  orthogonal  field  remains  excluded  because  of  the  low  dif¬ 
fusion  speed  of  that  field  relative  to  the  highly  conductive  plasma  mass. 
Thus,  the  high  conductivity  prevents  the  vacuum  field  from  entering  the 
plasma  when  the  curren!  .  ?  large  and  the  orthogonal  field  arising  from 
the  vacuum  field  of  the  accelerator  coil  geometry  cannot  be  effective.  The 
argument  is  different  for  a  single  coil  accelerator.  There  the  inefficiency 
is  due  to  the  conductivity  not  being  high  enough  to  limit  diffusion  of  the 
orthogonal  field  (which  also  drives  the  current  in  that  geometry). 
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SECTION  IV 


PULSED  BIAS  FIELD  EXPERIMENTS 

A  simple  argument  has  been  formulated  to  show  the  dependence  of 
efficiency,  for  a  self-field  accelerator,  on  the  diffusion  limited  force 
density.  The  current  density  j  at  any  point 


j  =  (T  (E  -  vB) 

where  a  is  plasma  conductivity,  E  is  electric  field  parallel  to  j,  vis  plasma 
velocity  and  B  is  magnetic  field  perpendicular  to  the  direction  of  motion. 


Thus 


V  B 


and  the  input  power  density  at  the  local  point  under  consideration  is 


.2 

jE=J5-+  (jB)v 


(6) 


It  is  equal  to  the  sum  of  the  local  ohmic  heating  and  the  product  of  the 
j  B  force  times  the  plasma  velocity.  The  latter  term  is  the  rate  of  work 
done  on  the  plasma.  This  consists  of  the  work  of  deformation  plus  the 
increase  of  kinetic  energy.  At  one  extreme,  the  "slug  model"  regards 
the  accelerating  mass  as  rigid;  at  the  other  extreme  the  "snowplow 
model"  assumes  impact  acceleration  in  which  the  work  of  deformation  is 
equal  to  the  kinetic  energy  gained.  The  ratio  of  the  two  terms  on  the  right 
side  of  Equation  (6)  is  simply  the  magnetic  Reynolds  number  for  the  accel¬ 
erator.  This  is  seen  by  writing  the  current  density  in  terms  of  the  Curl 
of  H, 


1  B 


(7) 


where  6  is  the  current  layer  thickness.  Then  from  Equations  (6)  and  (7) 
the  ratio  of  the  work  term  to  the  dissipation  is 


vBff  vBffp6 
j  "  B 


fi(r6  V 
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This  is  the  ratio  of  the  mass  velocity  to  the  magnetic  field  diffusion 
speed,  which  is  the  magnetic  Reynolds  number.  When  this  ratio  is  of 
the  order  of  1  to  3,  from  50  to  25%  of  the  power  goes  into  ohmic  dissipa¬ 
tion. 

In  the  bias  field  concept,  the  magnetic  field  in  Equation  (6)  is  the 
bias  field  which  is  perpendicular  to  the  acceleration  direction  z, 
while  that  in  Equation  (7)  is  parallel  to  the  acceleration  direction. 

The  ratio  of  the  work  term  to  the  dissipation  term  is  then 

z  z 

The  effective  Reynolds  number  is  then  increased  by  the  ratio  of  the  bias 
field  to  the  field  whose  diffusion  drives  the  plasma  current. 

4.  1  BIAS  FIELD  WITH  TWO-DIMENSIONAL  CURRENT  ACCELERATOR 

The  bias  field  concept  was  tested  by  adding  a  pair  of  bias  field 
coils  to  the  two-dimensional  current  accelerator.  The  field  coil  geometry 
is  shown  in  Figure  67.  Five  concentric  coils  were  wound  in  grooves 
machined  in  nylon  blocks.  They  were  then  encapsulated  in  urethane  and 
covered  with  a  5-mil  sheet  of  Mylar,  As  shown  in  Figure  68,  these  bias 
coil  blocks  were  placed  against  the  active  faces  of  the  two  large  accelera¬ 
tor  coils.  Before  construction  of  these  coils,  it  was  verified  that  a  good 
plasma  discharge  could  be  produced  in  spite  of  the  greatly  increased 
separation  between  the  current  drive  coils  and  the  plasma.  Figure  69 
shows  the  measured  bias  field  values  for  a  range  of  capacitor  voltages. 

The  field  and  current  configurations  for  this  experiment  are 
shown  in  Figure  70.  The  parallel  accelerator  current  sheets  produce 
the  antiparallel  fields  across  the  plasma  current  sheet;  the  gradient  of 
this  field  is  essentially  the  current  density.  The  bias  field  acting  to  the 
left  interacts  with  the  plasma  current  to  produce  the  upward  acceleration 
force. 

Figure  71  shows  both  the  circuit  current  and  the  plasma  current 
density  for  a  10-kv  accelerator  capacitor  voltage  and  a  static  fill  density 
of  350  millitorr  of  argon.  A  bias  field  circuit  capacitance  of  10  pF  was 
used  with  a  5i2  series -connected  damping  resistor  to  prevent  ringing. 
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Figure  67.  Bias  Coil  Geometry 


Figure  68.  Photograph  of  Accelerator  Coils 


Figure  69.  Measured  Bias  Field  for  Different  Voltages 
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Figure  70.  Current  and  Field  Configuration  in  Accelerator 
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Figure  71  Circuit  Current  and  Plasma  Current  Density  vs  Time 
at  Midplane  (350  m  torr  argon,  sweep  5  psec/cm) 

For  a  bias  capacitor  voltage  of  10  kv,  the  stored  energy  was  500  joules 
and  for  the  same  voltage  the  accelerator  capacitors  stored  900  joules. 
Using  the  bias  field  of  0.67  Webers/M  from  Figure  69  and  the  plasma 
current  density  from  Figure  71,  the  force  density  near  the  center  of  the 
bias  coil  and  midway  between  the  accelerator  coils  varies  with  time  as 
shown  in  Figure  72.  The  integrated  impulse  density  was  142  Ns/M  .  If 
the  volume  over  which  this  impulse  density  exists  was  lim.ited  to  the  inner 
diameter  of  the  bias  coil  and  a  typical  current  layer  thickness  of  1  cm 
were  used,  then 


Al  =  0.  069  NS 


Figure  72.  Time  Variation  of  Force  Density  at  Midplane 
(350  m  torr  argon) 
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and  with  the  conservative  assumption  that  represents  the  energy  stored 
in  both  capacitors, 


n  _  g  M  _  /O.  Q69\ 

I  2  J  \2M  1400  ' 
sp  o 


2.4  X  10 


/sec 


If  the  impulse  estimated  in  this  simplified  manner  were  applied  to  a  mass 
large  enough  to  limit  the  1^^  to  2  500  seconds,  the  efficiency  would  be  60%. 

4.2  COAXIAL  ACCELERATORS 


The  coaxial  geometry  is  naturally  suited  to  a  practical  thruster. 
Figure  7  3  shows  the  plasma  current  drive  coil  concept.  As  in  all  of  the 
previous  two  coil  accelerators,  the  plasma  current  is  driven  by  diffusion 
of  antiparallel  magnetic  fields  into  the  current  sheet.  In  the  coaxial  geo¬ 
metry  these  fields  are  developed  outside  of  the  inner  coil  and  inside  of 
the  outer  coil.  The  antiparallel  fields  are  not  present  in  the  vacuum 
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Figure  73.  Perspective  Showing  Field  Annihilation  Cur¬ 
rent  Drive  Coils 
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fields  of  these  coils  but  result  from  the  combination  of  those  fields  with 
the  fields  associated  with  the  plasma  current.  Thus  diamagnetic  current 
sheets  must  form  close  to  coil  surfaces  early  in  the  discharge  for  the 
coaxial  coils  to  be  an  effective  current  generator.  The  development  of 
such  a  discharge  is  one  of  the  practical  problems  of  the  coaxial  accelera¬ 
tor. 

Another  practical  problem  is  to  shape  the  magnetic  fields  in  the 
plasma  volume  such  that  the  average  mass  motion  will  be  axial  only. 

There  are  two  parts  to  this  problem.  One  is  that  the  normal  0-pinch 
force  of  the  outer  coil  be  balanced  by  the  inverse  0-pinch  force  of  the 
inner  coil  so  as  to  prevent  radial  drift  of  the  plasma.  The  other  is  that 
the  bias  field  be  essentially  radial. 

Two  axially  symmetric,  bias  field  accelerators  were  built  and 
tested,  both  using  coaxial  current  drive  coils  and  an  additional  set  of 
coaxial  coils  to  produce  the  radial  bias  field.  One  of  them  used  solenoidal 
coils  with  a  steadily  increasing  wire  spacing  to  produce  a  nearly  uniform 
radial  fringing  field  in  the  plasma  region,  while  the  other  used  four  coils 
arranged  to  produce  a  "cusp”  bias  field.  The  bias  circuits  were  much 
slower  than  the  drive  circuits  in  both  cases  so  that  the  bias  field  was 
fully  developed  before  firing  the  drive  circuit.  To  prevent  large  voltages 
being  induced  in  the  bias  coils  by  the  fast  transients  in  the  drive  coils, 
shorted  turn  shields  were  placed  close  to  the  bias  coils. 

The  first  axially  symmetric,  bias-field  accelerator,  designated  as 
Mark  I,  had  a  plasma  volume  1  cm  wide  by  6  cm  long  with  a  midradius 
of  5  cm.  The  coil  arrangement  is  shown  in  Figure  74.  The  current  drive 
coils  are  uniformly  spaced  but  a  constantly  increasing  spacing  is  used  for 
the  bias  field  coils;  the  spacing  ratio  of  adjacent  turns  is  a  constant  such 
that  the  spacing  of  the  last  pair  of  wires  is  ten  times  that  of  the  first  pair. 
The  plasma  volume  is  offset  as  shown  to  use  a  region  of  the  fringing  fields 
of  the  bias  coil  where  the  radial  component  is  reasonably  uniform  and  the 
axial  component  is  small  compared  to  the  radial  one. 

The  individual  coils  are  shown  in  the  exploded  view  of  Figure  75. 

The  inner  shield  can  be  seen  but  the  outer  one  is  not  because  it  was  used 
to  support  the  outer  bias  coil  during  winding.  The  assembled  accelerator 
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Exploded  View  Showing  Bias  and  Drive  Coils 
and  Shield  Assembly  Positions 


Figure  7  5 
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is  shown  in  Figure  76.  'Fhe  second  accelerator,  Mark  II  is  shown  in 
Figure  77.  A  photograph  of  the  accelerator  firing  is  shown  in  Figure  78. 
Shown  in  these  photographs  is  a  cross-shaped  arrangement  of  four  thin 
sheets  of  epoxy-fiberglass  that  serve  to  interrupt  the  plasma  current  as 
the  plasma  leaves  the  region  of  uniform  radial  bias  field.  Rogowsky 
probe  measurements  showed  that  this  technique  worked  well.  The  annular 
space  occupied  by  plasma  is  2  cm  wide  by  3  cm  long  with  a  mean  radius 
of  4  cm.  The  inner  and  outer  bias  coil  radii  are  2,  5  and  5.  5  cm  and  their 
axial  separation  is  6  cm.  The  plasma  volume  is  centered  axially  and 
radially  between  the  bias  coils. 

The  long,  narrow  slot  of  the  Mark  I  accelerator  was  found  to  be 
impractical  because  of  severe  probe  interference  with  the  current  sheet. 
Also,  in  this  accelerator  the  inner  and  outer  drive  coils  had  the  same 
number  of  turns  (one  turn  each).  This  may  have  resulted  in  a  radially 
inward  motion  of  the  plasma  during  acceleration  because  of  the  greater 
inward  "pinch"  force  from  the  outer  drive  coil  as  compared  to  the  outer 
"inverse  pinch"  force  from  the  inner  drive  coil.  These  difficulties  were 
overcome  in  the  Mark  II  design  by  using  the  wider  and  shorter  slot  and 
by  using  two  turns  for  the  inner  coil  with  one  turn,  as  before,  for  the 
outer  coil.  Probe  measurements  showed  the  axial  magnetic  field  to  be 
zero,  at  the  mid- radius,  at  the  time  of  maximum  current. 

A  peculiarity  of  the  discharge  was  observed  during  probe  measure¬ 
ments  with  the  Mark  II  accelerator  which  appears  to  indicate  that  an 
inductively  applied  bias  magnetic  field  is  not  feasible.  Figure  79  shows 
oscilloscope  traces  of  the  radial  magnetic  field  at  5  cm  radius  for  several 
axial  positions.  The  circuit  current  (lower  trace)  is  also  shown.  Zero 
lines  for  both  traces  were  recorded  for  convenience.  The  initial  bias 
field  of  0.  29  Webers /M^  can  be  seen  before  the  discharge  starts.  It  is 
seen  that  the  bias  field  is  reduced  nearly  to  zero  at  the  time  of  maximum 
accelerator  current,  near  the  closed  end  of  the  accelerator,  and  it  is 
increased  to  more  than  twice  the  initial  bias  value  at  the  accelerator  exit. 
This  effect  is  caused  by  the  current  pulse  induced  in  the  shorted  turn 
shields  when  the  accelerator  is  fired.  This  current  pulse  prevents  the 
magnetic  fields  from  the  accelerator  drive  coils  from  coupling  into  the 
bias  field  ceils,  as  intended.  However  these  closed  loop  currents  are 
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Figure  77.  Photograph  of  Mark  II  Accelerator 


Figure  78.  Photograph  of  Mark  II  Accelerator  Firing 
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Figure  79 

Radial  Magnetic  Field  at  5-cm 
Radius  in  Mark  II  Accelerator 
for  Different  Axial  Positions  _ 
(Initial  Bias  Field  =  0.29  W/M  , 
Imax  =  3900  amps,  2  psec/  cm) 


85 


strong  enough  and  close  enough  to  the  plasma  volume  to  have  a  first  order 
effect  on  the  radial  magnetic  field  while  the  plasma  current  is  running. 

Figure  80  shows  the  axial  distribution  of  the  first  half-cycle  maxima 
of  current  density,  bias  magnetic  field  and  force  density.  The  static 
fill  density  was  1/2  torr  of  argon;  a  2  pF  capacitor  charged  to  1  5  kv 
was  used  in  the  accelerator  circuit  (225  joules).  This  figure  shows  that 
the  shield  currents  not  only  had  a  large  effect  on  the  radial  magnetic 
field,  but  also  strongly  affected  the  plasma  current  as  well.  Thus  the 
plasma  current  was  large  where  the  field  was  low  and  low  (even  reversed 
at  the  exit)  where  the  field  was  large. 

The  reason  for  this  coupling  effect  can  be  seen  by  referring  to  the 
schematic  of  the  currents  and  magnetic  fields  of  the  Mark  II  accelerator 
in  Figure  81  .  The  drive  coils  produce  the  plasma  current  shown  and  the 
bias  coils  produce  the  outward  directed  magnetic  field  B^.  The  pulsed 
currents  in  the  shorted  turn  shields  have  the  polarity  shown,  i.  e.  ,  they 
have  the  same  direction  as  the  plasma  current.  The  incremental  radial 
fields  caused  by  these  currents,  AB^  in  the  figure,  oppose  the  bias  field 
at  the  closed  end  and  aid  it  at  the  open  end.  This  effect  is  readily 
apparent  when  it  is  recognized  that  all  currents  induced  by  the  accelerator 
coil  current  must  be  in  the  opposite  direction  as  the  coil  current.  Another 
equally  apparent  observation  is  that  each  pair  of  concentric  shield  coils 
acts  like  a  flat  coil  accelerator  which  drives  the  plasma  away  from  it  by 
inducing  a  plasma  cu”rent  counter  to  its  own  current.  Thus  the  transient 
fields  induced  by  the  shield  currents  produce  an  electric  field  opposite 
to  that  induced  by  the  drive  coil  electric  field,  thereby  reducing  the  plasma 
current  density. 

The  Mark  I  accelerator  configuration,  shown  in  Figure  82,  is  also 
affected  adversely  by  its  shield  currents  although  the  details  are  slightly 
different.  Again,  the  shield  currents  are  in  the  same  direction  as  the 
plasma  current  and  therefore  induce  an  opposing  electric  field.  The 
AB^  effect  is  opposite  to  that  in  the  Mark  II  accelerator.  The  radial  field 
is  increased  near  the  closed  end  and  reduced  near  the  open  end. 
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Figure  82.  Schematic  of  Currents  and  Magnetic  Fields 
in  Fringing  Bias  Field  Accelerator  -  Mark  1 


SECTION  V 


CONCLUSIONS 

1  .  For  the  plasma  accelerators  studied  in  this  program,  having  an 
argon  fill  density  of  0.  5  torr  and  a  magnetic  Reynolds  number  of 
the  order  of  3,  acceleration  by  the  direct  conversion  of  magnetic 
field  energy  to  plasma  kinetic  energy  was  not  observed. 

2.  Electron  temperatures  as  high  as  5  eV  were  produced  for  short 

4 

periods,  with  corresponding  plasma  conductivities  of  2x10  per 
12M.  Good  agreement  between  measured  conductivity  and  classical 
theory  (scattering  by  Coulomb  collisions  only)  was  found. 

3.  The  current  sheet  collision  developed  an  18-fold  increase  in  ion 
density  that  bounced  back  to  launch  magnetoacoustic  waves  back 
toward  the  drive  coils.  Just  outside  of  the  neutral  line,  at  80%  of 
the  coil  radius,  a  strong  density  gradient  developed  due  to  differen¬ 
tial  arrival  times  that  produced  a  large  velocity  early  in  the  collision. 

4.  Acceleration  by  self-fields  alone  was  inefficient  in  all  of  the  devices 
studied.  This  was  shown  to  be  an  unavoidable  consequence  of  the 
low  magnetic  Reynolds  number  of  the  small  scale  accelerators  of 
interest  for  electric  propulsion. 

5.  Efficient  acceleration  was  demonstrated  for  a  configuration  having 
orthogonal  bias  field  and  drive  field  coils. 

6.  Attempts  to  adapt  a  pulsed  bias  field  to  a  coaxial  accelerator  geo¬ 
metry  were  not  successful  because  of  currents  induced  by  the  drive 
coils  in  the  shield  rings  that  were  required  to  decouple  the  coaxial 
bias  and  drive  coils. 

7.  Some  new  technique  is  needed  to  supply  a  strong  bias  magnetic  field 
before  high  efficiency  plasma  acceleration  can  be  achieved  in  a 
practical  thruster. 
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